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ABSTRACT 


Motivated by the desire to put a free electron laser (FEL) weapon on a ship, 
the FEL and the related process of Compton backscattering are studied. The theme of 
the majority of this work is the interaction of the Gaussian optical mode with a beam of 


relativistic electrons. 


Classical FEL theory is reviewed in Chapter II. Simulations based on the 
classical theory are used in Chapter Ill to study a proposed 1 kW (kilowatt) infrared 
FEL. In Chapter IV, simulation is used to study the problem of electron beam/optical 
mode overlap in an ultraviolet (UV) FEL. A new concept, the FEL with a short 
Rayleigh length, is studied in Chapter V. The idea is tested on the UV FEL, then used 


to design and simulate a megawatt-class FEL for ship self-defense. 


An analytical calculation of the Compton backscattering of laser light is 
performed in Chapter VI. A quantum electrodynamics (QED) formalism is used to find 
the spectrum and angular distribution of photons scattered out of a Gaussian optical 


mode by relativistic electrons. 
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l. INTRODUCTION 


Today's US Navy is the most powerful the world has ever known. It can 
position aircraft, ships, and submarines within striking distance and project its power 
deep within the borders of almost any country. Yet today’s ships are arguably more 
vulnerable to attack from mines and anti-ship cruise missiles (ASCMs) than ever 
before, as a result of the Navy’s shift from blue to brown water. Ship defense against 
ASCMs is mainly a problem of reaction time; that is, there is not much time (usually 
seconds, depending on the missile’s speed and the ship’s detection range) between 
detection of the missile and impact of the missile with the ship. In this case a self- 
defense weapon which travels at the speed of light would be most useful. With 
continuing improvements in ASCM speed, stealth, and maneuverability, the problem 
gets even more difficult for existing defenses. A directed energy weapon (DEW) could 


be the solution. 


Experiments in missile destruction and atmospheric propagation have shown 
that a megawatt (MW) class free electron laser (FEL) is the directed energy weapon of 
choice for ship self defense (SSD). Megawatts are needed in order to deposit the 
required energy on the target in a short amount of time. More power probably would 
not shorten this time because thermal blooming, which results from atmospheric 
absorption, will degrade propagation. The FEL is needed because of its inherent 
tunability and designability, characteristics no other laser system can boast. The laser 
weapon's wavelength must operate in very narrow atmospheric absorption windows in 
order to propagate the beam to the target, avoiding problems like thermal blooming. 
The FEL has another advantage over other laser systems in that its gain medium is a 


vaCuum. 


While SSD is seen as the mission area where laser weapons can get their 
Start, a MW-class FEL aboard a ship could be an awesome offensive weapon as well. 


It was demonstrated in the 1970s that it is easier for the laser weapon to destroy 


crossing targets than incoming targets. It is possible that the precision and power of 
the FEL could also be used against aircraft and their pilots, satellites, ballistic missiles, 
other ships, and even personnel. Unlike the proposed Arsenal Ship, which merely 
packs a ship with more missiles, a shipboard directed energy weapon would truly bring 
about a revolution in military affairs, which would hasten the obsolescence of the 
missile. In this work we study the FEL and related phenomena, contributing to a field 


of scientific research with the goal of putting a DEW on a ship. 


Several topics are presented on the interaction of laser beams with relativistic 
electrons; this includes the free electron laser. In the next chapter, an overview of 
classical FEL theory necessary for understanding the ensuing chapters is provided. 
We describe the components of the FEL and then find the emission wavelength using 
basic physical arguments. Using the Lorentz force and Maxwell’s equations, we derive 
a set of two equations called the FEL pendulum and wave equations which describe 
the FEL interaction. Using these equations, we find expressions and describe the 
mechanisms for FEL gain. Modifications of the FEL interaction due to short optical 
pulses or the Gaussian optical mode are considered in the last two sections of the 


chapter. 


Original work simulating a proposed one kilowatt (1kW) upgrade to the 
FIREFLY FEL at Stanford University is presented in Chapter Ill. This upgrade would 
represent a one hundred-fold increase in FEL average power over that which currently 
exists. This work was originally presented at the 16th International FEL Conference 
and published as “Simulations of the Stanford FIREFLY 1kW free electron laser” in 
Nuclear Instruments and Methods in Physics Research (NIM) A 358, ABS 44, (1995). 


In Chapter IV is presented original research conceming optical mode 
transverse effects in ultraviolet free electron lasers (UVFEL). Specifically, it is the 
UVFEL proposed by the Laser Processing Consortium (LPC), a group based at the 
Thomas Jefferson National Accelerator Facility (TJNAF), that is simulated. The results 


were presented at the 17th Intemational FEL Conference and published as 


"Transverse effects in UV FELs" in NIM A 375, ABS 61, (1996). 


New and original work on a new concept in FEL design, the short Rayleigh 
length FEL, is presented in Chapter V. The concept is first tested on the LPC UVFEL. 
Those results were presented at the 18th International FEL Conference and will be 
published in 1997. We then use this concept to design and simulate a MW-class FEL 
oscillator suitable for ship self-defense. This design and the simulation results were 
presented at a workshop at TJNAF entitled "Navy MW-Class SSD FEL Concepts” 
which was sponsored by the Navy High Energy Laser office at the Space and Naval 


Warfare Command (SPAWAR). 


In Chapter VI we present new and original calculations of Compton 
backscattering of laser light from electrons. This process has been proposed as a 
means of generating x-rays for many applications, including medical imaging, and j- 
rays for a so-called y-y collider. Our calculations are the first attempt to incorporate the 


finite size of the laser beam into the usual scattering problem. 





ll. CLASSICAL FREE ELECTRON LASER THEORY 


A. THE FREE ELECTRON LASER 


The Free Electron Laser is a device which transforms the kinetic energy of a 
beam of relativistic electrons into coherent radiation (Ref. 1). A transverse, static 
magnetic field causes the electrons to wiggle and thereby radiate electromagnetic 
energy (Ref. 1). The radiation can either be stored in an optical cavity and amplified 
over many passes, as in the FEL oscillator, or amplified in a single pass, as in the FEL 
amplifier (Ref. 1). In this work, we will be concerned only with the oscillator 


configuration. 


The major components of the FEL oscillator are the electron accelerator, 
undulator, and optics. The accelerator delivers a beam of relativistic electrons to the 
FEL with energies anywhere from a few MeV to GeV depending on the wavelength of 
laser light desired. There are many different types of accelerators suitable for use in 
an FEL, but we will only be concerned with the radio-frequency linear accelerator (RF 
LINAC). The RF LINAC accelerates pulses of electrons, called micropulses, using 
either standing or travelling waves of an electromagnetic field (Ref. 2). One 
micropulse can be accelerated within each Serenata of the RF field, called a 
“bucket.” For example, if the frequency of the RF radiation is f = 500 MHz, buckets 
are separated by a distance c/f =60 cm. Electron micropulses are created and 
inserted into the accelerator by the injector (or electron gun). In our example, if the 
injector also operated at 500 MHz, it would fill every RF bucket and micropulses would 
be separated by 60 cm. The current in the micropulse is called the peak current I, 
and the average current is calculated using the duty factor (the ratio of the micropulse 
length to the distance between buckets). If we take / = 500 A and a micropulse length 
I, = 0.6 mm (/,/c = 2 ps), then the duty factor in our example accelerator is 0.001. 


The average current is then 7 = 500x0.001 A = 0.5 A. 


The undulator is a set of permanent magnets or electromagnets arranged to 
provide a periodic, transverse, magnetic field of wavelength A, along which the 
electrons traverse. The polarization of the undulator’s magnetic field can be either 
helical or linear. The world’s first FEL used a helical undulator, but the linear 


undulator is more common now (Ref. 3). Undulators can range in length from 1m to 


25m. 


The mirrors of the optical cavity (in the oscillator) store the light which is then 
amplified over successive passes of the electrons through the undulator. Total optical 
cavity losses, including the outcoupled light, are quantified by the quality factor Q, so 
that losses per round trip of the optical pulse are (100/Q)% (Ref. 3). The optical cavity 


is generally centered on and several times longer than the undulator. 


B. THE RESONANCE CONDITION 


With a basic knowledge of the physics of the FEL interaction we can define 


some important parameters and find the wavelength of the emitted radiation. 


Depicted in Figure 1 is the FEL interaction: the interaction of relativistic 
electrons with the undulator magnetic field, B.. and the radiated electric, E,, and 
magnetic, B,, fields. At the top is a schematic of an FEL showing the electron beam 
path through the undulator, the orientation of the undulator magnets, and the mirrors. 
The middie sketch is a blow-up of one wavelength A, of the undulator. The undulator 
magnetic field B oscillates in the yz-plane, causing the electrons to wiggle in the xz- 
plane. At the bottom of the figure is a blow-up of one half of one wavelength A of the 
radiation field. A portion of the electron trajectory is drawn at an angle to the z-axis. It 
is the transverse component of the electron velocity, that which projects on the xy- 
plane, which couples to the radiation electric field, allowing energy to be transferred 


between them. 
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Relativistic electrons interacting with an 


undulator magnetic field, B,, and radiated electric. E,, 


and magnetic, &,, fields (after Ref. 4). 


In order for the radiation in the oscillator cavity to grow, the electrons must lose 
energy to the light wave on each pass through the undulator. The Lorentz force 
shows that this energy loss is dy/dt « B E. where B = V/c, V is the particle velocity 
and the relativistic factor y = (1 — B*)"'? (Ref. 5). Since E, is transverse, that is, in a 
plane perpendicular to the z-axis, only a transverse component of B couples to the 
cavity radiation. Driven by the undulator magnetic field, that transverse component is 
of the form B, « Re (Ke"°7) where K is a constant called the “undulator parameter," 
Ky = 2m/Ay, and A, is the undulator wavelength. The electric field is of the form 
E, « Re (E,e!\*2-°'9)) where k = 2n/A, 4 is the radiation wavelength, w = kc, and 6 is 
an optical phase factor. Therefore, we find that dy/dt « Re (KE, tho ~ ON 6 io), It is 


evident from these preliminary considerations that the electron “~ponderomotive" phase 


defined as 

C=(k + ko)z - ot (1) 
is fundamental to the FEL interaction. For relativistic FELs, this phase represents the 
position of the electron within a radiation wavelength. 


The time it takes the electrons to traverse the undulator is L/cB, = 10°. 
Rather than work with time scales on the order of 10° seconds, we introduce a 


dimensionless time t, defined as 


—> 


2 cB t C 


a 
—_— 


Sen an 





(2) 


where L is the undulator length and cB, =v,=c in the relativistic limit. As the 
electrons traverse the undulator t ranges from 0 to 1. We can now define the electron 


phase velocity as 
v(t) = ot) = Lik ee ee (3) 


The phase velocity measures the rate of change of the electron's phase relationship 


(C) to the optical field. 


In order that the rate of transfer of energy from the electrons to the radiation 
field is maximized, the wiggling electron and the electric field must be nearly resonant. 
Looking at Figure 1, we can explain this as follows: when the electric field is pointed 
in the x-direction, the electron trajectory needs to have a component in the x-direction. 
When the electric field phase increases by n and is pointed in the (-x)-direction, the 
electron trajectory must also have a component in the (-x)-direction. This will be the 
case if exactly one wavelength of light passes over the electron every time the 
electron traverses exactly one wavelength of the undulator. This is the "resonance 
condition.” lf one wavelength of light passes over the electron for every undulator 
period the electron traverses, the electron will maintain the same phase relationship 


with the light. Resonance thus implies that the phase velocity v(t), must be zero. 


Invoking the resonance condition, we set v(t) = 0 in (3) and solve for A with the 
result ’ = Ao(1 — B,)/B,. From the definition of y, we find y* =1- B? — 82. In the 


helical undulator (as determined later), Be = K*/°. Solving for B, we get 


1+K° 
2y 


for y >> 1. Substituting (4) into the expression for A above we find 


Bz =1- — (4) 


1+ K- 


ey 


which is the wavelength of light emitted by the electrons in the FEL. This expression 





= ho (5) 


shows that we can change the wavelength by changing the electron beam energy 


ymc? or the undulator characteristics 4, or K. 


In the preceeding discussion we defined and then discussed resonance for a 
single electron. in a real electron beam, there are about 10'° electrons in the 
micropulse; all of them could not possibly be exactly resonant. Although we define 
resonance as v = C = 0, where ( ‘ ) represents differentiation with respect to t, we will 
find that in the low gain FEL (section E) it is desirable to operate at just above 


resonance v(0) = vo = 2.6. We can understand this by considering that about 10° 


electrons are uniformly spread over the length of an optical wavelength. If they are 
resonant, they maintain the same phase relationship with the light throughout the 
interaction even though the light is passing over them. Thus at resonance the 
radiation from the i” electron at C; will exactly cancel that from an electron at a phase 
(C; — =). This cancellation will occur for each electron. To avoid this cancellation, for 
net radiation, the electrons must radiate together, or "bunch" on the scale of the 
optical wavelength. This requires that ee and t >0. We will discuss this 


microscopic bunching in part three of the next section. 


C. ELECTRON DYNAMICS 


Electrons in the undulator interact with both the undulator field and the radiation 
field which has been set up in the cavity. In subsection 1, we find the electron 
trajectories in the undulator using the Lorentz force equation of motion. We then use 
these trajectories in the second subsection to deduce measures of electron beam 
quality. In the third subsection we assume that electron radiation has set up a 
classical radiation field in the optical cavity and then allow it to interact with the 


electrons. This will result in the pendulum equation, which describes bunching. 


Ae Electron Trajectories in the Undulator 


A linear undulator polarized in the y-direction causes electrons to wiggle in the 
xz-plane. In addition, since electron beams have a small radius and angular spread 
(which depends on the quality of the beam from the accelerator), electrons off axis in 
the y-direction are in stronger fields and are therefore focussed back toward the axis. 


This results in so-called betatron oscillations. 


The magnetic field of a linear undulator with no external focussing in the x- 


direction is given by 


B =[B,, B, B,) = B[0, sin(kpz)cosh(koy), cos(kpz)sinh(kgy)] (6) 


10 


where By is the amplitude of the undulator field (Ref. 3). Consider an idealized 
electron which is injected into the undulator perfectly (i.e., x(0) = B,(0) = 0) in the x- 
direction. The electron motion is governed by the relativistic Lorentz force equation of 
motion in the undulator field (Ref. 5): 

& = -eB xB (7) 
where e = !e| is the magnitude of the electron charge, and the electron’s momentum 
sao. = ymcRB, m is the electron mass, and c is the speed of light in vacuum. Note that 
cgs units are used. Substituting (6) into (7), noting that dy/dt = 0 since the magnetic 


field can do no work on the electron, we get the set of equations 


B 
B, = EBzsin(koz (t))eosh(koy(t)) - Bycos(koz(t))sinh(key(t))] (Bal 
eBy , 
By = eB xcos(koz(t))sinh(koy(t) (8b) 
Deh 
Bz = ra on Moma COS oF AU) (8c) 


where ( )=d/dt. The right side of (8a) is an exact time derivative and can be 


integrated immediately with the result: 
By = “2 cos(kyz)oosh(koy) (9a) 


where K = eBA,/2nmc® and B = Bo/V2 is the root-mean-square (rms) magnetic field 
strength. The dimensionless parameter K is the undulator parameter and is an 
important quantity in FEL physics, as shown in (5). The constant of integration has 
been set to zero indicating our assumption of perfect injection in this direction. The 
time dependence of the dynamical variables has been suppressed for brevity. 


Substituting (9a) into (8b) and (8c) we get 


CKaka 
By = a? (KoZ)sinh(2koy) (9b) 


11 


Chelan 9 
i sin(2kyZ)cosh*(Kgy) . (9c) 


~ 


Equations (9) are coupled, nonlinear and too difficult to solve analytically. In 





order to investigate the xz-plane motion we will now add the restriction of perfect 


injection in the y-direction also. That leaves us with 


V2K 





Be = Simic Oz) (10a) 

B,(0) = y(0) = 0 (10b) 
Chak, 

i 9 sin(2k oz). (10c) 


Now in the FEL, K=1 and y> 1, so that K/y « 1 (Ref. 3). This allows us to expand 
equations (10) for the trajectories, keeping only terms to order (K/y)*. To zeroth order 
in (K/y), z(t) =cB,t where cB, we take to be the average electron velocity in an 
undulator period. Using (4), we find cB, = c(1 - (1+K2)/2y° + ...)=c. Substituting 


z(t) = cB,t into (10c) to find the next higher order contribution we get 


Pen ee aeemer cane eae (11) 


Ako 


which is already to order (K/y)*. To find x(t), we will substitute the first term of (11) 


into (10a) since using the second would lead to higher order terms in x(f). We get 
If 
x(k = _»2K cin(kocB,t) a (12) 
Koy 


where in both (11) and (12) we have set B, = 1 in the denominator. 


We can rewrite (11) and (12) in terms of the dimensionless time t (See (2)) as 


x(t) =- 





K? 
ae sin(2nNt)+ ---> (13a) 
V2ry 


2 





z(t)=B,lt- - sO cant) t+ occ: (13b) 


V2 


where N = L/A, is the number of undulator periods. As expected, the x-component of 
the trajectory is the wiggling motion caused by the undulator with one oscillation per 
undulator period. The z-component contains the average velocity term and a term 
with a fast oscillation (twice per undulator period). The oscillation is superimposed on 
the average motion because as an electron with constant speed makes an excursion 
to the side, the z-component of its velocity must slow, leading to an oscillation in z- 
trajectory. An electron makes this excursion twice per linear undulator period, which 


causes the z-trajectory to oscillate twice per undulator period. 


Now that we have found the trajectories in the x- and z-directions (13), we will 
find the motion in the y-direction. The y-trajectories are called betatron oscillations 
and are caused by the focussing action of the undulator magnetic field. We need to 
use a more realistic beam in this direction, but maintain a small enough beam such 
that Kpy«1 and sinh(2kpy) = 2koy. If we then re-write (9b) in terms of t, the result is 

Be 2K*L°k§ 


(t) = 2 


Since a typical betatron oscillation takes several undulator periods to complete, we can 


cos*(2nNt)y(t) . (14) 


average (14) over several periods (Ref. 3). We are then left with 


“an aa ye _ Ve(y) 
y 2 ay 


which is the equation of motion for a simple harmonic oscillator. The solution, and 





(15) 


thus the y-trajectory, is sinusoidal. The harmonic potential is Vs(y) = wsy*/2 where 
Wy = KLko'y = 2nNK/y is the dimensionless betatron frequency. In a typical undulator 
N = 100, K = 1, y= 100, and we get w, = 2x, or one betatron oscillation along the 


undulator. 


We have only treated betatron motion in one direction, but there can be 
betatron focussing in the other direction as well. This could be provided by extemal 


focussing magnets, or by machining parabolic pole-faces on the undulator magnets 
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(Ref. 3). Neither the amount of focussing nor the electron beam quality have to be 


the same in both directions. 


2. Electron Beam Quality and Emittance 


By inspection of (15) and knowledge of Hamiltonian mechanics, we know the 
harmonic betatron oscillations have a constant of motion 
Hg = [py2(t) + y2(t)/2 = [wpy2(0) + y2(0)V2. It can also be shown that when an 
electron is injected off axis (y(0) = yo) or with an initial angle (y (0) = L 6, where 6, is 
the angle from the z-axis), its z-velocity is decreased by AB, = —Hp/L* relative to a 
perfectly injected electron (Ref. 3). Taking differentials of (3) we find Av = LkAB, since 


k > ky. Therefore, an imperfectly injected electron experiences a decrease in phase 


velocity given by 


kH 

B enn oe 2 

Ay = -—— =- K“k 0 16 
V L are oYo + ¥8y) (16) 





where we have used (5). 


In a real electron beam from an RF LINAC there can be =10'° electrons per 
micropulse, each of which has a phase velocity change away from the average given 
by (16). Numerically we don’t simulate all 10'° electrons, but instead distribute many 
electrons in both injection position and angle. They are distributed as a Gaussian, 
which has the distribution function 


f(q) = exp(-g*/20°) (1%) 


V2n 6 


where o is the standard deviation and q is the vanable being distnbuted (Ref. 3). 


Using (16) as a guide we define the standard deviation of the angular distnbution as 


awe 2nN/o* (18) 


re 1+ K? 


where 8 is the root mean square (rms) initial angle from the z-axis. The standard 
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deviation of the position distribution, also called the dimensionless electron beam 


radius, is 


Oona = Fo (19) 
where fp, is the rms initial distance from the z-axis (Ref. 3). The significance of the 
factor VAL/r, the characteristic optical mode size, will be discussed in section G of this 
chapter. From (16) we see that 


_ 2nNK?kérG 


is the standard deviation of an exponential spread in phase velocity caused by the 
initial Gaussian spread in position characterized by (19). The Gaussian spread in 
injection angles characterized by (18) also leads to an exponential distribution in 
phase velocity with the same standard deviation (Ref. 3). The exponential distribution 
arises because the phase velocity of an electron decreases no matter what the 
direction of the (Gaussian distributed) initial position or angle. The distribution function 


for the exponential is given by 


bia rs q<0 (21) 


where q is the variable being distributed, i.e., v;, the phase velocity of the i” electron 


(Ref. 3). 


The injection angles and positions of the electrons in the beam are governed 
by a quantity called the “emittance,” which is a measure of the quality of the electron 
beam. The rms emittance is defined aS t~n. = Mp8, and is a fixed quantity throughout 
the entire length of the undulator (Ref. 2). That is, one cannot change the radius of 
the electron beam without affecting its angular spread and vice versa. Solving (18) 


and (19) for 8 and 7, and substituting these into the definition of emittance yields 


fans = VT hE (22) 
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where we have defined the dimensionless emittance € as 
2Sonan - (23) 


This dimensionless emittance will be used extensively in Chapters IV and V. 


o The Pendulum Equation 


We now consider the interaction of the electrons with both the undulator and 
radiation fields. We have already found the macroscopic electron trajectories in the 
undulator. Our goal in this subsection is to find an expression which describes the 
evolution of bunching. We will use the field of a helical undulator in this section as it 
makes the problem mathematically easier. At the end we will state a prescription for 


transforming the pendulum equation to that for a linear undulator. 


The magnetic field on the axis of a helical undulator is given by 
B, = By[cos(koz), sin(koz), 0] (24) 
and the radiation field present in the cavity can be represented by 
E = E(x,y,z,t)[cosw, —siny, 0] (25a) 
B = E(x,y,z,t)[siny, cosy, 0] (25b) 
where E(x, y, Z, ft) is the wave amplitude and w = kz — wf + o (Ref. 3). 


To describe this problem of the interaction of electrons with the undulator and 


radiation fields, there are five coupled equations: 


= -e(E + Bx B) (26) 
oy. Oe 
dt mae E a 


1 


1 Fare 
T= . 


but only four unknowns: B,, By, Bz, and y, so we can neglect the z-component of (26). 


(28) 
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The fields E and B in (26) and (27) are the sum of the radiation and undulator fields. 


Substituting (24) and equations (25) into the x- and y-components of (26), we get 
d e . . 
<(1B,) = --2{E(1 ~ B,)[cosy, -siny, 0] + B,Bol-sin(koz), cos(koz), O}] (29) 


where we have used Pp = ymcB, and E = E(x,y,z,t). The right side of (29) contains 
two terms, the first of which is proportional to (1 — B,) and is surely small for relativistic 
electrons. We will neglect this term and keep only the dominant undulator term (i.e., 
that which is proportional to By). We can then integrate (29) and set the constant of 


integration to zero indicating perfect injection. The result is 
K : 
By = es coS\ 07: sin(KoZ), 0] (30) 
where K = e€ByAp/2nmc* is the undulator parameter since B= Bp in the helical 


undulator. 


As mentioned earlier, E is transverse and couples with only the transverse 


components of B, that is, (30). Therefore substituting (80) into (27) we are left with 


(_ EK 
Y= Tit cos(¢ + 9) (31) 


where € is given by (1) and 9 is the phase of the optical field as discussed above. 


As we noted earlier, (28) can be arranged as y*=1-B*-f%. If we 


substitute (30) for B? we get 


1— 82 
tee on 





which shows that a change in y does not result in a change in #,, only in B,. Taking 


the time denvative of both sides of (32) results in 


1. ¥BeB, 


(33) 
Y 1+ K? 


Construction of y'Y from (31) and equating that with (33) results in 


Me 


3s 14K? (K+ko)C eEK 
C= 2 er cos(f + 4) (34) 








where we have used B, = C(k + Ko)c. We will now make approximations using 


k > Ky and B, = 1, and substitute for k using (5) to get 


_ 2k oeEK 


C Pr imac ; (35) 


We can make (35) dimensionless by setting d?/dt® = (c/L)*d?/dt°, which leaves us 
with 


oo 


C = lalcos(C + 6) (36) 


which is the pendulum equation. We have defined the dimensionless optical field 


amplitude lal as 


al= 


~ me? (37) 


The pendulum equation appropriate for the linear undulator can be found by the 
substitution K — K[Jo(&) — J1(&)] in (37) (Ref. 3). Here J,(&) and J,(E) are Bessel 
functions of the first kind of order O and 1 and — = K/2(1 + K°). 


Equation (36) shows that electrons with phase —7/2 < (C + 6) < w2 have a >0 
while those with w2<(€+06)< 37/2 have a <0. This difference in phase 
acceleration shows how electrons which are initially distributed randomly from 
(C + o) = -w2 to ($ + o) = 32/2 are driven together, i.e., bunched, near (f + >) = 2 by 
the FEL interaction. In the next section, we will see what affect these bunched 


electrons have on the optical field. 


D. THE FEL OPTICAL WAVE EQUATION 


Knowing through the pendulum equation that the electrons will be bunched by 
the FEL interaction, we now find what affect those bunched electrons have on the 


classical radiation field. We use Maxwell's equations driven by the electron current to 


18 


determine the form of the FEL optical wave equation (Ref. 3). 


Electrons in the undulator emit spontaneous radiation because they are made 
to accelerate. After approximately one pass, there is enough radiation in the optical 
cavity to take on a classical form (Ref. 3). Over the next several passes this 
spontaneous emission spectrum is narrowed. In many cases, the final FEL emission 
spectrum, or linewidth, is determined by the length of its optical pulses through the 
Fourier transform (Ref. 3). Short optical pulses result in broader spectra than long 
optical pulses. At this point in the operation of the FEL, the radiation wavelength is 
given by (5). Since we know that the emission spectrum of the FEL is narrow (it is a 
coherent light source), it must be that the amplitude and phase of the radiation field in 
the cavity are slowly varying in both time and space. We will invoke the slowly-varying 
amplitude and phase approximation in the following derivation of the FEL optical wave 


equation (Ref. 3). 
The fields in the cavity of an FEL with a helical undulator are given by 


equations (25). Since the amplitude and phase are slowly varying in time, we can 


write the vector potential for these fields as 


ius suze ive (38) 


where € = [-/, 1, 0] is the polarization vector and y = kz - wf +. The potential for 
the physical fields is given by the real part of (38). The inhomogeneous wave 


equation for A is given by 





V2 - 43a og, (39) 


where J is the driving current density (Ref. 5). The slowly-varying amplitude and 
phase approximation can be written as inequalities: E « kE, 9’ « ko, E <wE and 
 < wd with ( )’ = o&()/dz (Ref. 3). In the wave equation, this approximation has the 
effect of favoring first derivatives over all other terms. Substituting (38) into (39) and 


Carrying out the 2- and t- derivatives we are left with 
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Ve te 2K 2 i E(x,y,z,t)e!%yzt) _— 288 Tages — wt) (40) 
Cc 








where V? = 0*/dx* + d°/dy” and é* is the complex conjugate of é. 


The driving current J is given by 
J (x,y,z) =-ec ¥ (x — x) &y - yj) &(z - z) By =-ec Y SPF) B, (41) 
f f 
where 06(x) is the Dirac delta function, 7 is the position vector of the field point 


(x, y, Z), F is the position vector of the i” electron, and the index / runs over all 


electrons in the micropulse. Substituting (30) and (41) into (40) and dividing by 2ik 


results in 


Ee’? = perci® 





(42) 





i xr 0 1d 
L203 det ot 


where we have used (1) and E = E(x,y,zZ,). 





We made the left side of (42) slowly varying by invoking the slowly varying 
amplitude and phase approximation. The right side is not slowly varying: it contains 
many Fourier components due to the presence of the delta functions which are each 
infinitely broadband. This means that while the left side describes an optical field that 
varies slowly, the right side can describe an optical field which varies arbitrarily fast. 
To smooth the right side we must replace the sum over all electrons by an average 
over some sample electrons. Since the fields in the cavity are slowly varying over 
several wavelengths by assumption, we only need to sample the electrons in one 
section of the micropulse one wavelength of light long and use the local charge 
density to weight the average. The result is 

| | Ay? + 


ae i 
| 4x oz cat 





a 12 





Ee'® = -2neK p(z - cB,t) < © oma as ep,» (43) 


where p(z - cf,t) is the charge density of the traveling electron pulse, and <---> is 


an average over sampled electrons. 
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We can write (43) in a much simpler form if we change variables to Z = z —- ct 
and use the dimensionless time t. The coordinate Zz then follows a point on the 
optical pulse traveling at speed c from which the electron pulse slips back at speed 
c(1-—8,). In this coordinate system (let Z > Z) we find that L(d/dz + 1/cd/dt) = d/dt. 
We can make (43) dimensionless by making the replacement x,y — Va/AL x,y where 


Vi.L/n is the characteristic optical mode size. The result is 


a cee 
-iv+3 





-it 
Ee'® = -2neKLp < = > (44) 


where we have assumed a long electron pulse so that p can be considered constant. 


Also, if the FEL efficiency is small so that changes in y are small, (44) takes the form 


-iV} + ala =—j<e!S (45) 





where a = lale’®, lal is given by (37), and the dimensionless current density j is 


ae soe | (46) 


To get the wave equation for an FEL with a linear undulator, we make the same 


substitution as before, K —- K(Jo(E) — J4(E)), in both (37) and (46) (Ref. 3). 


Now we will consider the physics contained in (45). First, we see that the 
growth of the field amplitude d!a|/dt depends on both the current / ‘and the wavefront 
Curvature through the diffraction operator ve This unnecessarily complicates our 
discussion, so we will assume for the moment that the light is constructed of plane 


waves so that E(x,y.z,t) ~ E(z,t) and Via — 0. Then (45) takes the simple form 
pS icou al (47) 


We see that if s = 0, then a = 0 and the field does not grow. Also, if there is no 
bunching, the electrons are randomly distributed throughout each wavelength of light, 
<exp(-/f)> = O and again the field does not grow. Equating the real and imaginary 


parts of (47) leads to 


eu 


dd lea: 
sae j<cos(f + o)> (48a) 


ce 
la| aie j<sin(C + o)> (48b) 


which shows that the field amplitude |a! grows when the electrons are bunched near 


(¢ + >) = m whereas the phase 9 is driven when bunching is near (C + 6) = 7/2. 


The wave equation (45) and the pendulum equation (36), when coupled, allow 
the field and the electron phase space positions to evolve self-consistently. The 
pendulum equation describes the bunching of electrons in the presence of the 
undulator and radiation fields, while the wave equation shows how much the radiating 
bunch changes the optical field. The two equations form a feedback loop which under 


the right starting conditions leads to gain, our next topic. 


E. GAIN 


Now that we know how the electrons affect the optical field and vice versa we 
can translate this interaction to gain. FELs can be separated into two regions: low 
currentlow gain defined by / « 1 and high current/high gain defined by / > 1 (Ref. 3). 
As it turns out, the low gain behavior we will analyze has a broader range of 
applicability, up to at least f/=1. An FEL with ;=100 will exhibit high gain 
characteristics. The difference between the low and high gain regions is in the growth 
rate of the optical field a (see (47)). We will find that the gain mechanisms in these 


two regions are different and therefore warrant a distinction between the two. 


Gain for either high or low current is defined as 


G= lal*- ag _ Energy Lost by Electrons (49) 
als Optical Field Energy 


where the last equality arises from conservation of energy. Here ao is the optical field 


strength seen by the electrons upon entering the undulator. Using (3), (4), and (5) it is 
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not difficult to show that the change in electron phase velocity due to a change in 


energy Ay/y or wavelength emitted AA/A is 


Ne ann (50a) 
or 
Ae anno* | (50b) 


so that the energy lost by one electron can be expressed as 


(Vo = Vv) 2 


Ae = mc 
nN y 


(51) 


where E = ymc? is the electron energy. The number of electrons in a volume of the 
optical mode dV is pFdV, where the filling factor F is the ratio of the cross sectional 
areas of the electron beam and optical mode. The optical energy contained in that 
same volume dV is E*dV/4n (Ref. 5). After some algebra we can write (49), letting 
j — jF, as 


o 


eV £L eyo-v> (52) 
d( 


which is valid for both low and high current. An average over the phase velocities v; 


in the sampled portion of the electron beam is indicated by < --: >. 


In the two subsections which follow we will discuss gain in the low current and 
high current FELs, respectively. We will find expressions for the single pass gain in 
both cases as a function of t. In FEL design, low gain usually implies an oscillator 
configuration while high gain (e.g. j/ = 10‘) implies an amplifier in which electrons 
execute only one pass before the light is outcoupled. In this case, the dimensionless 
power P = |a!* evolves as P = (1 + G)P, and a large G ensures P > Po. After n 
passes in the oscillator P, = (1+ G)"P,>) so although G is small, over many passes 


the light is amplified until P, > Po. 
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In the following analysis the wave equation is in the form (47). In section G we 


will discuss the effects of diffraction on the FEL interaction. 


le The Low Gain FEL 


In exact analogy with a simple pendulum, the tool used to analyze the slow 
electron dynamics (bunching) in the FEL is their (C,v) phase space, one 2z section of 


which is one wavelength of light long. 


In the low gain FEL, / <1 implies through (47) that changes in the field are 
small. Multiplying both sides of (36) by 4 and rearranging terms leads to 


Cale 
—}|— - lal = 53 
a | i alsin(t + ’ 0 (53) 
and we recall that v = C. This has a simple solution 
1/2 
v(C) = 2/2 + 2lalsin(¢ + 0) (54) 


where Hy = vé/2 - lalsin(Co + 0) is a constant of the motion for each electron starting 
at (Co, Vo) and determines which phase space path v(C) it will traverse. These paths 
are plotted in Figure 2. For electrons starting on resonance vo = 0, (54) predicts fixed 
points at (Co + 6) = 7/2 and 3n/2. Analysis shows that the fixed point at 7/2 is stable 
and surrounded nearby in phase space by harmonic orbits (Ref. 6). . The fixed point at 
3/2 is unstable and defines the path characterized by Hy = lal, called the separatrix, 


which is 

v, = tv2Ial i + sin($ + 6)}"2 . (55) 
The separatnx separates phase space into regions of closed and open orbits and has 
a peak to peak height of 4\/a! (Ref. 6). It is plotted in bold in Figure 2. 


In Figure 3 we have distributed 20 resonant (i.e., Vo = 0) electrons uniformly 
through phase space and show their evolution as t1:0—> 1. At t=0, the electron 


position is gray and steadily darkens as t > 1. The phase space simulation solves the 
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wave and pendulum equations numerically with / = 1 and initial field amplitude a, = 3. 






FEL PHASE SPACE PATHS 





_——————— 


-1/2 C 3/2 


Figure 2. The phase space paths v(C) traversed by the 


electrons. The separatrix is drawn bold. 


since the electrons are resonant, they start and remain inside the separatrix on closed 
orbits. One can clearly see in Figure 3 that the electrons become bunched at € = n/2, 
as predicted at the end of section C. Even though the electrons are bunched, 
however, gain is nearly zero because the electrons are bunched at the wrong phase 
for field amplitude growth. We can see in the figure that as many electrons move up 
in phase space and gain energy from the light wave (see (50a)) aS move down in 
phase space and lose energy to the light wave. Something must be done to upset this 


balance tn order to get a net gain. 


It is interesting to note that Figure 3 shows a very small amount of gain. This 


is due to the slight phase shift given by (48b) which causes the bunch to move just 


ZS 


slightly toward (¢ + 0) = 72. 






| *** FEL Phase Space Evolution *** 








Figure 3. Phase space simulation of low gain FEL with 


electrons starting on resonance. Gain is negligible. 


To get an analytic solution for the gain in a low current FEL, expand (36) in 


weak fields la! = ag « x with the result 


V(t) = Vo + ae 'sin(o + VoT) - sin(So) + (56) 
Vo 


2 
a 
: 4 [cos(25o + 2vot) - cos(2éo) + COS(VoT) — 1 — VotSIN(So)COS(So + vo “tl 
L 


where >) = !a(0)! and (0) = O (Ref. 3). When we perform the average in (52) only 
the term proportional to af in (56) survives and the result is 


2 - 2cos(vot) - a 57) 


60-1 a 
0 
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We see that the gain at the end of the undulator (t = 1) depends only on the initial 


phase velocity and current density when the system is in weak optical fields. 


In Figure 4 is the result of a simulation which solves the pendulum and wave 
equations numerically for many different initial phase velocities and plots the final gain 
G(t = 1) and phase shift o(t = 1) versus vo. The vg axis can be viewed as either a 
function of electron energy or emitted wavelength through equations (50). In the 
figure, Vo ranges from -12 to 12, /=1, and ag =0.001 <x. The gain spectrum is 
anti-symmetric in agreement with (57) and has a maximum of G = 0.13 at vo = 2.6. 
The phase curve shows that the shift is maximum for resonant electrons, as discussed 


above, and is quite small where the gain is large. 


Oe ren oe 









***k Gain and Phase Curves **** 
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Figure 4. Gain and phase spectra for low gain FEL. 


Gain is maximum at Vo = 2.6. 
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We have in Figure 5 another phase space simulation using the same 
parameters as in Figure 3 except that vo = 2.6 in order to maximize gain. In this case 
only about half of the electrons are trapped in the separatrix while the other half are 
situated on open orbits. A clear bunch of electrons can be seen at ¢ = x which affords 
maximum gain with a small phase shift ¢. Compared with Figure 3, more electrons 
lost energy to the light wave than absorbed energy from the light wave. It is also true 
that if the electrons started at vy = —2.6 the result would be maximum absorption as 
more electrons would be absorbing energy than losing energy. This behavior could 
have been predicted from studying Figure 2: since all electrons travel nearly the same 
distance along their orbit (from (36)), starting them higher in phase space (i.€., Vo > 0) 
results in electrons with G) < 7/2 traveling more or less parallel to the ¢ axis, while 


those with C, 2 7/2 travel down, losing energy. 






*** FEL Phase Space Evolution *** 
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Figure 5. Phase space simulation of low gain FEL with 


electrons starting at vo = 2.6 for nearly maximum gain of 


G = 0.12 with ao = 3. 
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2. The High Gain FEL 


In the high gain FEL, j > 1 and changes in the optical field are not small. 
Energy conservation cannot be used to determine the gain in this case because there 
is a substantial change in both the optical field amplitude and phase. It is still possible 
to obtain an analytic expression for the gain in weak optical fields; we merely state the 


result: 
G(x) = Sexpl( Sj)" (58) 


where vo is set to zero (Ref. 3). Unlike the low gain FEL (57), here the gain is 
exponential in t. We can ascertain the dependence on vy by using the same gain and 
phase simulation used in Figure 4. Shown in Figure 6 is that simulation with J = 100 


and a) = 0.001 « x. 


**** Gain and Phase Curves **** 


j=100 a.=0.003 
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Figure 6. Gain and phase spectra for high gain FEL. 
Gain is maximum at vo = 1.4, but is still substantial at 


Vo = O where the phase shift o = w2. 
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The maximum gain is G = 71 at vo = 1.4, closer to resonance and much higher than in 
the low gain case (Figure 4). Note that the phase shift at maximum gain is quite large. 
The gain spectrum is nearly symmetric, and substantial gain exists at resonance 


where the phase shift = 7/2. 


To understand the mechanism for high gain, a phase space simulation using 
J =100, ag =1, and vyg=0 is shown in Figure 7. The electrons which start at 
resonance drive the optical phase o, represented in phase space by the shifting 
Separatrix, and attempt to bunch at ¢= 7/2 as discussed earlier when considering 


Figure 3. 


*** FEL Phase Space Evolution *** 





Figure 7. Phase space simulation of high gain FEL with 
electrons starting at resonance. Final gain G = 62 = e“ 
and the phase shift 6 = 7’2. The separatnx is shown 
growing with the field strength !a! and shifting left as the 


phase 6 grows. 
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It is the phase shift @ = 2/2 which allows gain as the relative electron bunch phase 
(¢ + 6) = x drives the optical field amplitude. This leads to exponential gain (linear on 


the logarithmic plot of Figure 7) as predicted by (58). 


To summarize, when / < 1 we must start the electrons at vg > O in order that 
they bunch near (¢ +) = and gain results. For an FEL with / > 1, electrons at 
Vg = O drive the optical phase with the high current until the bunch is situated near 
(C+) =7. In the case of low gain, the optical phase shift is unimportant and small; in 


the case of high gain it is a crucial aspect of the interaction. 


3. Gain Degradation Due to Electron Beam Quality 


In section C.2 of this chapter, we discussed transverse electron beam quality 
and deduced two measures, 6g and o,,., from the emittance. The longitudinal beam 


quality has a measure 


Og = 4nN es (59) 


which is the standard deviation of a Gaussian spread in energy (see (17)), and Aymc* 


is the rms energy difference from the resonant energy (Ref. 3). 


Consider the low gain spectrum in Figure 3 again. The full width of the main 
positive gain portion at half the maximum is approximately m and the distance from 
maximum gain to maximum absorption is approximately 2x. Therefore, if the spread in 
phase velocities of the electron beam Av=n7, we should expect serious gain 
degradation as many electrons will be injected into the undulator with phase velocity 
corresponding to less than maximum or even to negative gain. Inspection of (50a) 


and (16) shows that Av has contributions from og, 6,,, and Gg. 


In Figure 8 is a phase space simulation using the same parameters as those 
used in Figure 5 except that a small energy spread og = 1 has been included. Only 


the final phase space positions of the electrons have been plotted. The electrons are 
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spread through phase space and bunching has been degraded. Gain has been 


reduced to G = 0.10. 






*** FEL Phase Space Evolution *** 
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Figure 8. Phase space simulation of low gain FEL with 
energy spread 6g =1. Gain is reduced to G =0.10. 


Compare with Figure 5. 


In the high gain spectrum of Figure 6 we note that the gain spectrum bandwidth 
for j = 100 ts approximately 2x. In general, this width is proportional to j'® in the high 
gain region (Ref. 3). Since there are no regions of absorption and the spectrum ts 
about twice as wide as that of the low gain example, we would expect the high gain 
FEL to be more resistant to beam quality effects. However, when increasing the 
undulator length, beam quality measures og. (59), dy. (18), and o,,. (20) increase 
proportional to N. This means that the quality of the electron beam ts magnified by a 


long undulator. Since ; ~« IN°, where / is the peak current, high gain is most easily 
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and often obtained by using a long undulator. In addition, higher peak current / from 
the accelerator leads to poorer beam quality, so that either way the designer tries to 
maximize gain the problem with beam quality gets worse. When trying to get all one’s 
light amplification in a single pass, little gain degradation can be tolerated. Indeed this 
may become one of the major drawbacks of using an FEL amplifier as a high power 


weapon. 


In Figure 9 is a phase space simulation using the same parameters as those 
used in Figure 7 except in this case Gg = 1. Only the final separatrix and electron 
phase space positions have been plotted. The electrons are spread throughout phase 
space and bunching has been degraded. Gain has been reduced from G = 62 to 


G10: 
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Figure 9. Phase space simulation of high gain FEL with 
energy spread 6g =1. Gain is reduced to G = 50. 


Compare with Figure 7, where G = 62. 
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In both Figures 8 and 9 we have only included the energy spread og 
contribution to beam quality. Clearly the degradation would be even worse if we also 


included the effects of emittance. 


FE: SHORT PULSES 


lt was mentioned during our derivation of the wave equation that the electrons 
slip back relative to the light pulse at speed c(1-—8,). From the definition of 
resonance, we know that the distance between the leading edges of the electron and 
light pulses at t=1 is NA. This distance is called the "slippage length.” If the 
electron pulse length is on the order of the slippage length, the pulse is considered 
short (Ref. 3). In this case we cannot make the long pulse assumption in (44) 
because each section of the light pulse sees the charge density changing as the 


electron pulse slips past. 


In an optical resonator of length S, it takes time 2S/c for the optical pulse to 
make one round trip. Let the desynchronism d be the distance between the light and 
electron pulses at t=O on each pass (Ref. 3). We normalize d and the electron 
pulse length o, = /,/(NA) to the slippage length. At exact synchronism, unbunched 
electrons slip back past the light pulse as t : 0 — 1, bunching occurs, but only the rear 
of the light pulse gets amplified. On successive passes the delay in bunching 
becomes longer and the amplification moves further back on the light pulse. Since 
only the rear of the light pulse is being amplified while the losses (as determined by 
the resonator quality factor Q) are the same at all points in the pulse, its centroid 
moves at a speed slower than c. This is called “lethargy,” and ultimately will lead the 
pulse out of the region of gain (that is, it won't overlap the electron pulse) and the 


pulse amplitude will vanish (Ref. 3). 


To compensate for the effects of lethargy a small amount of desynchronism 


must be added. This is accomplished by reducing the path length S that the light 


pulse must travel by adjusting the length of the optical cavity. In this case, 


d =—-AS/(NA) (Ref. 3). 
Simulations of an FEL using picosecond electron pulses will be discussed in 


Chapter lll. There we will explore the effects of desynchronism on the power 


evolution. 


G. FEL WITH GAUSSIAN OPTICAL MODE 


In section E on gain we discussed the interaction of the electron beam with 
light approximated by plane waves. An optical cavity bounded by spherical mirrors 
supports Hermite-Gaussian (and other) modes, which are solutions of (45) instead of 
(47) and accurately describe the laser beam inside and outside the cavity (Ref. 7). 
The optical beam can be constructed from any number of these modes, but the one 
most desired in a weapons laser for its propagation characteristics is called the 


fundamental Gaussian mode, or just Gaussian mode. 


The usual Gaussian mode as given in Ref. 7 can be written in terms of our FEL 


dimensionless parameters. The result is 


2g -riw?z,_ ida(t 
e 0 G(T) 








ag = 60 
Ciara (60) 
where 
aT ee 
$g(t) = -tan7“! | ———| + anmenauia (61) 
Zo Z6 + (T- Ty) 
is the phase shift of the Gaussian mode relative to that of a plane wave, 
2 
t-T 
w*(t) = 1 | = (62) 
Zo 


is the square of the mode radius, which is normalized to the mode waist radius wé. 


and 
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Zo=Wwé6 (63) 
is the dimensionless Rayleigh length (Ref. 3). The mode waist position is t,. Radial 
dimensions are normalized to VLA/n, the characteristic mode size. We see in (62) that 
Zo is the distance over which the mode area doubles. In a symmetric optical cavity 
with spherical mirrors a distance S apart and radii of curvature R, the dimensionless 
Rayleigh length is given by z)=(2SR -— S*)”/2L where L is again the undulator 
length (Ref. 7). The Rayleigh length Z) and the waist position +, completely 


characterize the Gaussian mode in the optical cavity. 


The amplitude of the Gaussian mode la! is a maximum at t=7, where 
w*(t)=1. It decreases away from the waist in both directions and is a minimum at 
both t= 0 and t = 1 (if t, = 0.5). At any time t, the amplitude decreases in the radial 
direction as a Gaussian, hence the name. The phase of the Gaussian mode $< (61) 
contains two terms. The first term provides a phase shift, called the Guoy phase shift, 
relative to the plane wave, and switches sign at the waist (Ref. 7). The second term in 
(61) decreases the phase shift for positions off-axis which creates the wavefront 


curvature necessary for focussing the beam at the waist. 


In the electron phase space, the term a)/w(t) causes the height of the 
separatrix to increase as t:0— 17, and then decrease as t: ty > 1. The effect is 
proportional to w~“(t) (from our knowledge of the separatrix height) and complicates 
the phase space paths. The phase of the Gaussian mode (61) is continuously shifted 
as t:0-— 1. The shift has the sign opposite to that caused by the FEL interaction in 
(48b). It is therefore represented in phase space by the separatrix shifting to the right 
as T:0—> 1. Looking at Figure 2 again, we can predict that the electrons will have to 
Start at a phase velocity vy. > 2.6 in order to compensate for the shift and still bunch at 
the correct phase for gain. (Compare this to the case of the high gain FEL where the 


phase velocity had to be lowered in order to bunch at the correct phase.) 
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In order to observe the effects of the Gaussian mode on the FEL interaction, 
assume Zp > 1, nearly plane waves, and expand (61) and (62) to first order in (1/Z,). 
The phase becomes 6g = -(t - Ty)/Zp, and w(t) = 1. The pendulum equation (36) in 
the presence of the Gaussian mode is t = laglcos(C + dg) which for Z,) > 1 becomes 


vig -r7/Z Tw 1 
C = aoe cos|Cy + —— + |Vo - — 
Zo Zo 


(64) 





when expanded in weak fields to first order in a, (Ref. 3). Comparing (64) with (36) 
we see three distinct differences. First, the amplitude of the phase acceleration a 
decreases for electrons off axis because they see a weaker field strength 
lal + ayexp(-r?/z,). Second, there is a shift in the initial phase C, of all electrons in 
the amount 1,,/Z>. This is inconsequential since the electrons are initially distributed 
uniformly throughout phase space. The last difference is in the phase velocity, which 
is decreased by Av=1/Zp. This leads to a wavelength change of AA/A = 1/(2nNZp) 
from (50b). This also requires an increase in the initial phase velocity for maximum 
gain to vf = 2.6 + 1/Z,) in the low gain FEL. This effect is caused by the optical 
phase shift (61) and is shown in Figure 10. The figure is the result of a simulation 
similar to that in Figures 4 and 6 except that the finite size of the optical mode is 
included. Here the Rayleigh length Z)>=5 and maximum gain is found at 


Vo = 2.9 = (2.6 + 1/5). 


The presence of the Gaussian optical mode produces measurable effects in the 
pertormance of the FEL. This is the theme of the bulk of this work, including the final 
chapter on Compton backscattering where we calculate the photon emission rate for 


electrons in the Gaussian field of a laser. 
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_*** FEL 3d simulation, single-pass gain *** 
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Figure 10. Gain spectrum for low gain FEL with Rayleigh 
length Z)>=5. Maximum gain is achieved at vo = 2.9, 


shifted from v, = 2.6 by the Gaussian mode. 
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ll. SIMULATIONS OF FIREFLY 1KW FEL 


A. BACKGROUND 


The FIREFLY (Far Infrared FEL) 1kW FEL at Stanford University was proposed 
in 1994 as an upgrade to the existing Infrared FEL at Stanford University (Ref. 8). 
With four new 10 MeV accelerator structures and a 5 MeV electron injector, an 
electron beam with energy 45 MeV could be generated. The electron injector 
operating at 11.8 MHz would produce an average current of /=1 ma, and the 
average power in the electron beam would be P, = 1V = 45 kW. A good measure of 
the maximum efficiency of an FEL, |.e., how much energy one can extract, is the gain 
spectrum bandwidth. Considering Figure 4 and (50b), we find that Ay/y = 1/(2N). If 
we assume that we can extract the full Ay/y = 1/(2N) from the beam with N = 25, then 
the average power out of the cavity would be P = P,/(2N) =900 W. This average 
power would be almost 100 times greater than the world’s most powerful FEL, and it 


was estimated this could be done for only one million dollars with existing technology. 


The proposed FEL and two variants with lower electron beam energies of 25 
MeV and 10 MeV were simulated. All other parameters were kept the same and are 


listed in Table 1. 


B. SIMULATIONS AND RESULTS 


Because all three FELs simulated have emission wavelengths A in the far 
infrared, the characteristic mode size VLWr is large enough that diffraction (i.e., the 
Gaussian mode) should not dominate the interaction. The slippage length NA, 
compared to the electron pulse length, is also large, so the effects of short optical 
pulses must be addressed. (Some basic consideration of short pulse effects was done 
in section F of Chapter Il.) To study short pulse effects, simulations following 


longitudinal multimodes were used. These simulations divide the optical pulse into a 
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number of fixed longitudinal field sites, or bins, a(z) (Ref. 3). The simulations are 
called "longitudinal multimode" because dividing the optical pulse into longitudinal 
sites a — a(z) is equivalent to a — a(k,) in Fourier space; the number of sites 2 is 


equivalent to a number of modes k,. 


Beam Energy 


Peak Current 

Micropulse Length 

Undulator Type Linear Linear Linear 

K 1.1 1.1 eat 

N 25 25 25 

Xo 6.0 cm 

r 156.7 um 
7.2 
0.2 





Table 1. Parameters for FIREFLY 1kW FEL (column |) 


and two variants (columns II and Ill) used in simulations. 


The slower moving electrons slip back past the optical pulse to sites 2 — t and interact 
with a number of field sites corresponding to the slippage distance NA (recall that Z is 
normalized to the slippage distance). Conversely, each field site a(z) interacts with a 
length N2. along the electron micropulse as the micropulse slips back. We generalize 


the pendulum (36) and wave (47) equations to include the slippage: 
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00 


G zat = la, Icos(o 74+ + o;) (65) 


an [eee : (66) 


As tT: 0 — 1 in increments of At in the simulation, (65) and (66) show how the field 
changes at site z are caused by the electron dynamics at the micropulse site z + s At, 
where s = 0, 1, 2, ...,N; is an integer and N, = 1/At is the total number of time steps, 


and vice versa. 


Contained in Figure 11 its the output from a simulation of the 1 kW machine, 
whose parameters are listed in column | of Table 1. The dimensionless current 
density is / = 1.7, the micropulse length is o, = 3.2, cavity losses are determined by 
@= 100, and the desynchronism d =0.03 is large (see Chapter II section F for 
definition of dc). 


| **** FEL Pulse Evolution **** 


j=1.7 O,=3.2 Q=100 d=0.03 
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Figure 11. Longitudinal multimode simulation of 


FIREFLY. 
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The parabolic electron pulse is drawn in the lower left window at t= 0 and at t= 11 to 
illustrate slippage. Above the electron pulse is plotted the optical field amplitude 
la(z,n)l at each pass n. The scale ranges from zero in black to the maximum in 
white; a contour line is drawn at the median field amplitude. Above the field evolution 
is plotted the final optical field amplitude !a(z,n = 2000)! in black. At the bottom 
center is drawn the gain spectrum (57). Above this is the optical power spectrum 
evolution P(v,n), with the final power spectrum P(v,n = 2000) drawn on top in black. 
The power spectrum P(v,n) is the Fourier transform of the square of the optical field 
la(z,n)l? at each pass n. The phase velocity v in the power spectrum is used in the 
sense of (50b); that is, it represents the wavelength A so that P(v,n) shows the 
linewidth of the radiation. The total power evolution P(n) is plotted at lower right in 
the figure. Above this is plotted the electron phase velocity distribution f(v,n). The 
distribution f(v,n) is a projection of the final electron phase space positions (C,v) onto 
the v-axis at each pass n. As such, v in this plot is to be understood in the sense of 
(50a); that is, it represents the electron energy ymc* so that f (v,n) is interpreted as an 


energy spread. 


Having described all the plots in Figure 11, we now interpret their features. 
The large desynchronism has caused the optical pulse !a(z,n)! to be artificially long 
by advancing the front edge of the pulse ahead of the electrons. The electron pulse 
amplifies only the trailing edge of the optical pulse on each pass. This leads to 
reduced coupling; the peak field strength is only |a! = 25. As evidenced by the power 
spectrum P(v,n), the long, smooth optical pulse leads to a narrow linewidth. The 
electron phase velocity distribution f(v,n) does not evolve any further after the first 
few hundred passes. At saturation, the electrons reach the bottom of the separatrix in 
phase space (Ref. 3). This means that the electron distribution f(v,n) should have a 
width of about 4v/al = 18 at saturation, the peak-to-peak height of the separatrix at 
la! = 20, the final field amplitude. In the figure, the width of f(v,n) is about 18, 


showing good agreement. The power P(n) at each pass rn is obtained by integrating 
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over the optical pulse la(z,n)l*. The power evolution P(n) plotted at lower right 


shows that a steady state is reached after approximately nm = 700 passes where 


Plotted in Figure 12 are the results of several simulations like that above. The 


Steady state dimensionless power is plotted versus the desynchronism. 
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0 0.005 0.01 0.015 0.02 0.025 0.03 
d 
Figure 12. Plot of outcoupled dimensionless power 


versus desynchronism for FIREFLY. 


The maximum on the curve at d = 0.0005 corresponds to an actual output average 
power P = 900 W., the desired output, but requires the injector to operate at 23.6 MHz, 
or double the rate discussed earlier. Evidently the short pulses lead to reduced 
coupling between the electrons and optical field, lowering FEL efficiency. Doubling the 
operating frequency of the electron injector was proposed in Ref. 8 as a follow-on to 


the 1 kW upgrade for the purpose of attaining 2 kW average power. Note that the 
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maximum in Figure 12 is rather sharp. Since the slightest change in d can produce a 
large change in the output, the FEL is less stable when operated in this region (Ref. 
=) 

The FIREFLY variant listed in column II of Table 1 has a shorter dimensionless 
electron pulse length o, than that in Figure 11. The output of a simulation with 
J = 2.9, 6, = 1, Q = 100, and d = 0.005 is shown in Figure 13. The features of the 


plots show the occurrence of the trapped-particle instability (Ref. 3). 
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Figure 13. Longitudinal multimode simulation for the 


variant of FIREFLY listed in column II of Table 1. 


This instability occurs when the field amplitude |a! becomes so strong that electrons in 
phase space traverse their entire path (see Figure 2) dunng a single pass through the 
undulator, executing a synchrotron oscillation with frequency v, = ‘la = 2x (Ref. 3). 


In real space, this oscillation represents motion superimposed on the electron 
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trajectory and leads to emission into sidebands vy + v, (Ref. 3). As these sidebands 
grow, the FEL power will increase and the fields will get stronger, causing more 
electrons to execute synchrotron oscillations, etc. At least four sidebands can be seen 
in the power spectrum P(v,n) in Figure 13. These sidebands modulate the optical 
pulse envelope !a(z,n)! as shown, and the increased field strength causes the wide 
electron distribution f(v,n). In this case, the instability has imposed a_ periodic 
structure on the optical field envelope |a(z,n)! and the power spectrum P(v,n) which 


results in oscillations in the total power P(n) over many passes n. 


Figure 14 plots the dimensionless power versus desynchronism from several 


simulations like that of Figure 13. 


900 
800 
700 
600 


500 





0 0005 OO OO01S 002 0025 0.03 
d 


Figure 14. Plot of dimensionless outcoupled power 
versus desynchronism for vanant of FIREFLY listed in 


column {Il of Table 1. 
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We see from the plot that Figure 12 with d = 0.005 lies in the transition region in the 
desynchronism curve, where the fields are strong enough to cause the trapped-particle 
instability, but where operation of the FEL is relatively stable. The peak in this curve 


at d = 0.003 corresponds to an outcoupled average power of P = 140 W. 


To further illustrate the effects of slippage, the example listed in column III 
represents FIREFLY operating at a much longer wavelength 2 = 156.7 um. The 
slippage length is therefore long and leads to a small value of o, = 0.2. Figure 15 
contains the result of a simulation using / = 7.2, 6, = 0.2, Q = 100, and the very small 


value of d = 0.0005. 






cee FEL Pulse Evolution **** 


j=7.2 O,=0.2 Q=100 qad=0.0005 







£ (V7 5) 





-4 2 =101l LY 101 = ee 


(2-0) SSSC«SY |P (a) 520 | 
tel | aaah ae | 
<< 


-4 z 4-101 Vv 101 0 n 6000 


Figure 15. Longitudinal multimode simulation for variant 


of FIREFLY listed in column Il! of Table 1. 


The small desynchronism means that, unlike in Figure 11, the light pulse .a(z,n)i has 


not been artificially lengthened. The optical pulse reaches steady state with a length 
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comparable to the electron pulse length shown in the plot at lower left. The fields are 
strong and the trapped-particle instability again leads to structure in the optical pulse. 
We should note here that the spikes to the left of the main optical pulse are only about 
one wavelength of light in length. Because this would violate the slowly-varying 
amplitude assumption in the wave equation (47), these may be meaningless. The 
short, modulated optical pulse results in a broad power spectrum P(v,n) (note the 
scale in v compared with that in Figure 13). The strong fields also result in a broad 
electron distribution. Maximum power P(n) = 900 is achieved at n = 1500 passes, but 
steady-state P(n) = 520 isn’t reached until n = 4000 passes which is predicted to take 
about 0.3 ms in FIREFLY. 


Figure 16 plots power versus desynchronism for column Ill of Table 1. 
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Figure 16. Plot of outcoupled dimensionless power 
versus desynchronism for varnant of FIREFLY listed in 


column Ill of Table 1. 
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We note that d = 0.0005, that used in Figure 15, is the location of the peak in this 


curve and corresponds to an average power of P = 17 W. 


C. CONCLUSIONS 


The operational characteristics of far infrared wavelength FELs driven by short, 
picosecond electron micropulses will be determined by short pulse effects. The 
desynchronism coupled with the cavity Q will affect the length and structure of the 


optical pulse. This in turn determines the width and structure of the power spectrum. 


Based on our simulations of the 1 KW FEL, as depicted in Figures 11 and 12, 
the goal of 1 kW average power can be achieved, but at double the electron injector 
duty cycle first proposed. More simulations should be done to substantiate this 
possibility. Doubling the injector duty cycle was proposed in Ref. 8, but as a follow-on 
to the 1 kW proposal as a means of achieving 2 kW. We also note that the trapped- 


particle instability will probably occur. 
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IV. TRANSVERSE EFFECTS IN UV FREE ELECTRON LASERS 


A. BACKGROUND 


In an ultraviolet free electron laser (UVFEL), the characteristic optical mode 
size VLA/n tends to be small because of the short optical wavelength A. For example, 
if A= 200 nm and L =3 m, VLA/x = 0.4 mm. This is also approximately the size of 
the electron beam. For the ultraviolet FEL, the effects of the spatial extent of the 


electron beam as compared to the optical mode are of concern. 


The UVFEL under consideration in this chapter is that proposed by the Laser 
Processing Consortium (LPC) at Thomas Jefferson National Accelerator Facility 
(TUJNAF) in Newport News, Virginia (Ref. 9). The LPC is a group of industry and 
university scientists and engineers interested in an efficient, reliable, coherent UV light 
source for processing materials. On the order of 100 kW of UV light is needed to 
process commercial products ranging from nylon to steel. The parameters for the 


UVFEL proposed in Ref. 9 are listed in Table 2. 


Beam Energy 200 MeV 
Peak Current | 270 A 
Micropulse Length 0.25 ps 
Undulator Type Linear 
K 0.95 
72 


N 
ho 
nN 





Table 2. Parameters for the LPC proposed UVFEL. 
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In this chapter we analyze the LPC UVFEL using transverse multimode 
simulations. These simulations are described in section B. We use these simulations 
in section C to study the effects of the optical mode size on FEL gain. In section D we 


study the effects of electron beam focussing on FEL gain. 


SB: TRANSVERSE MULTIMODE SIMULATIONS 


In Chapter Ill we extended the pendulum and wave equations in the z-direction 
in order to properly simulate longitudinal multimodes by dividing the optical pulse into 
z sites. To simulate transverse multimodes, we make the same transition in the x- 
and y-directions. The optical field is extended to sites in x and y so that a > a(x,y). 


Generalizing the wave equation (45) and the pendulum equation (36) gives (Ref. 3): 


‘ A : 
ae aE 2 a (x LY st) Sl exp(—/ O(X,YT))> x,y.) (67) 





oo 


C (xyt) = [alxy.n| cos|Soy.2) + oby.n)] (68) 


The transverse dimensions are normalized to the characteristic mode size YAL/x. The 


set of equations (67) and (68) are fully self-consistent in three dimensions (x,y,t). 


For the UVFEL represented in Table 2, we calculate an electron pulse length 
slightly greater than five slippage distances 6, =5.1. Therefore, it can be assumed 
that slippage will not be an important consideration, so that we may follow only one 


representative longitudinal site as in (67) and (68). 


In section G of Chapter Il, we discussed the fundamental Gaussian mode but 
mentioned that in general any number of Hermite-Gaussian transverse modes are 
present. In order to model the evolution of the FEL from weak fields through 
saturation, these higher-order modes are included as the FEL interaction amplifies and 
distorts the optical mode. Both the optical field amplitude and phase can be altered 


(Ref. 3). An example of this is shown in Figure 17. 
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Figure 17. Optical multimode evolution from weak fields 


through saturation. 
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Figure 17 shows the evolution of the FEL from weak fields ay = 0.0001 through 
saturation using / = 20 and a small energy spread og = 0.5. The Rayleigh length for 
the cavity is Z7=0.5. The top set of plots show the FEL oscillator after n = 12 
passes. The electron phase space plot shows no noticeable bunching and the final 
gain G=0.94 is in steady state (i.e, G(n) constant ) with power increasing 
exponentially. At right is shown the optical field amplitude between the mirrors 
la(x, t)l. The mode is clearly distorted from a fundamental Gaussian mode with the 
peak field (in black) located at t= 1 where single pass gain is the maximum. Moving 
down the set of figures, the electron phase space begins to exhibit bunching at 
no 24. The bunch descends along the v-axis, indicating energy loss, and near 
saturation at n = 28 it starts to ascend taking energy back from the optical field. At 
saturation (n = 40), the electron bunch is clearly moving up in phase space at t = 1. 
The motion of the electron bunch through phase space, caused by the strengthening 
fields, is attended by a drop in final gain G(t=1) as the electrons are absorbing light at 
the end of the undulator. This means that on each pass, G(t) peaks at t< 1 along 
the undulator and the peak field moves back along the undulator. Scanning down the 
page through the plots of la(x,t)l, we see the peak field strength moving back along 
the undulator. At saturation, the peak gain per pass is nearly at the middle of the 
undulator and the final gain is reduced to G = 0.0034. The optical field la(x, 1)! 


resembles the fundamental Gaussian mode at saturation where the gain is low. 


The evolution of the distorted optical mode to the Gaussian requires the use of 
the transverse multimode simulations described above. Multiple transverse modes are 
included in the simulation to construct the distorted optical fields of Figure 17 for the 
same reason it takes a large number of sine waves to synthesize a square pulse in 


Fourier analysis. 
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C. OPTICAL MODE SIZE EFFECT ON FEL GAIN 


It has been shown that the Rayleigh length which optimizes gain for a filament 
electron beam (i.€., 6g = 0) is Z) = 1/V12 = 0.3 (Refs. 3, 10). This is the Rayleigh 
length which maximizes the average filling factor F (averaged along the undulator 
length) by minimizing the volume occupied by the Gaussian optical mode around the 
electron beam. Reference 10 also discusses gain in a Gaussian mode for non-zero 
electron beam radius, but again using arguments based on the filling factor. The 
conclusion, not surprisingly, is that maximum gain is achieved when the Rayleigh 
length is chosen such that the electron beam fits inside the optical mode (recall 
Zy = wé). But in the case of the UVFEL, it may be necessary to have an electron 
beam bigger than the optical mode because of the small mode size VAL/n. In that 
case, one cannot quantify gain in terms of the filling factor, which may be greater than 
one. The three dimensional simulations described in the previous section are used to 


find the influence of the Rayleigh length on FEL gain. 


Figure 18 shows the result of a three-dimensional, transverse multimode 
simulation of the UVFEL parameters listed in Table 2. The corresponding 
dimensionless parameters are current density / = 30 with F is no longer included in j, 
Vo = 4.5 is the optimum initial phase velocity, the Rayleigh length is 2) = 0.3, 
A@p = 0.001 is the initial optical field amplitude, the cavity loss is Q = 5, o, = 0.45 is the 
dimensionless electron beam radius, dg, = 0.9 is the angular spread, the energy spread 
IS Gg = 1.3, and w,; = 1.0 is the betatron frequency. The dimensionless emittance 
¢ = 0.43 corresponds to an actual normalized emittance e, = 11x mm-mrad 
(Ey = ms). The plot at top left shows the evolution of the optical field amplitude at 
the end of the undulator with each pass |a(x,n)!. To the nght is plotted the final 
optical field amplitude at the nght-hand mirror la(x,y)!. The optical field between the 
mirrors |a(x,t)| at y = 0 is plotted in the center; the tick marks at t = 0 and t = 1 show 


the beginning and end of the undulator. Overlaid on this plot are the positions of a 
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random sample of the electrons used in the simulation. In this figure, o,/Wo = 0.8 and 


the electrons near the edges of the Gaussian mode see weaker optical fields and thus 


a degraded interaction (see (60)). 
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Figure 18. Transverse multimode simulation of UVFEL. 
Gain is degraded due to large electron beam 


radius/optical mode waist radius ratio 6,/Wop = 0.8. 


The plots along bottom show the electron phase velocity distribution at each pass 
f(v,n), the final phase space positions (¢,v) of a sample of electrons, and the power 
P(n) and gain In(1+G(n)) evolutions. The gain G(n) is plotted at the end of the 
undulator G(t=1) after each pass n. These plots show FEL gain in steady state (i.e.. 
G(n) approximately constant) with final gain G=0.75 and power increasing 
exponentially. Electrons in phase space are spread out due to 6g and oy So that 
bunching can barely be seen at ¢ = m even though the limit of weak fields !a! = n has 


been reached. 
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The results of many simulations like that in Figure 18 were used to plot gain 
normalized to the maximum one-dimensional gain, G )=0.135/F, where 


io 2 eel deeaicaine average filling factor, versus Rayleigh length (Refs. 3, 10). 


The result is plotted in Figure 19. 
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Figure 19. Plot of normalized gain versus Rayleigh 
length for UVFEL. Normalized gain is minimum where 


optical mode volume is minimum Zo = 0.3. 


The first feature to note in the figure is that for the UVFEL, 2, = 0.3 corresponds to the 
minimum normalized gain because the optical mode volume is minimized. The optical 
mode volume must be large enough to envelope the electron beam. Decreasing the 
Rayleigh length below 2, = 0.3 causes the mode to expand quickly away from the 
mode waist. This makes the mode large at the ends of the undulator, enveloping the 


electron beam, and increases normalized gain. Increasing the Rayleigh length trom 
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Zy = 0.3 makes the optical mode larger everywhere and increases gain until the mode 
gets so large, Z) = 0.9 in this case, that electrons see weaker fields everywhere and 
gain is reduced. For a larger electron beam, the optimum value of Z, for normalized 


gain would be larger than Zp, = 0.9. 


D. ELECTRON BEAM SIZE EFFECT ON FEL GAIN 


Now, given that the Rayleigh length of the optical cavity has been chosen, how 
is FEL gain affected by the choice of electron beam radius? Since ws, = 1.0 for this 
UVFEL, the electron’s betatron trajectories are nearly straight lines, allowing them to 
be focussed at the undulator center using external magnets. Keeping emittance 
€=6, V6, fixed, decreasing the electron beam radius at the waist increases the 
angular spread and degrades the interaction, so there must be a limit to how small 
one can focus the electron beam. In all simulations, the Rayleigh length is held fixed 


at a = 0.5: 


Figure 20 contains the output of a simulation of the UVFEL electron beam with 
€ = 0.43. In this example, we have reduced the electron beam radius to 6, = 0.3, 
which increases angular spread to 6, = 2.1. We also were required to increase the 
current density to / = 68 in order to keep peak current constant. The optimum initial 
phase velocity is vp = 4.3 and all other parameters are the same as those in Figure 
18. The phase space plot (¢,v) is dominated by the (exponential) angular spread 
which indicates 6, has been made too small. In the plot of |a(x,t)! in the center, one 
can see the electrons focussed to o, = 0.3 at the center of the undulator, but the large 
angular spread carries many electrons outside the optical mode at t= 0 and t= 1. 


This leads to gain degradation. 


Figure 21 plots the normalized gain versus electron beam radius for three 
different emittances: (a) ¢ = 0.28, (b) ¢ = 0.43, and (c) € = 0.6. (The simulation in 


Figure 20 is from curve (b).) We see in the curves that decreasing the electron beam 
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radius from o, = 1.0 increases gain as more electrons are enveloped in the optical 
mode. The beam radius can be made too small, however, causing gain to decrease 


because of the increasing angular spread. 


TRANSVERSE FEL MODEL 










2; ie ‘3 a Serbs 3 fet 
ee EA: ah = ode ea ates 
pes tof 


t 
Figure 20. Transverse multimode simulation of UVFEL. 
Electron beam focussed at t=0.5 but angular spread 


causes beam to expand outside optical mode at tT=0 


and 1. 


The optimum electron beam radius, and the peak gain at that radius, depend strongly 
on the beam emittance. With better emittance, curve (a) in Figure 21, the optimum 
electron beam radius 1s small, 6, = 0.3, and peak gain ts increased to G/Gp = 66%. 
Increasing emittance causes the optimum o, to grow larger and the maximum to 
become more broad until (curve (c)) the maximum gain ts only weakly dependent on 


G6, due to the large angular spread at any radius. This maximum gain is reduced to 


ay 


G/Gy = 10% at o, = 0.6. The expected emittance in the LPC UVFEL is e = 0.43, 
curve (b). We find that the optimum electron beam radius is co, = 0.4, slightly smaller 


than that proposed. 
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Figure 21. Plots of normalized gain versus electron beam 
radius for three different beam emittances: (a) ¢ = 0.28, 


(b) € = 0.43, and (c) e = 0.6. 


Ee: CONCLUSIONS 


The small characteristic mode size of the short wavelength UVFEL makes the 
electron beam/optical mode overlap a crucial consideration. Making the optical mode 
larger by either decreasing or increasing the Rayleigh length from Zp, = 0.3 results in 
increased normalized gain for a given electron beam radius. This is an important 
result for UVFELs. It shows that it is not necessary to have the electron beam 


enveloped by the optical mode throughout their interaction. With Rayleigh length fixed, 


08 


focussing the electron beam to a smaller radius will increase gain until the angular 
spread, which must increase as the electron beam radius is decreased, degrades the 
interaction and gain is reduced. With improvements in emittance, the optimum 
electron beam radius is smaller and the peak FEL steady-state gain is significantly 


increased. We find the optimum electron beam radius for the LPC UVFEL is o, = 0.4. 
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V. FELS WITH SHORT RAYLEIGH LENGTH 


A. MOTIVATION 


The problem of high power density on optical cavity mirrors is shared by both 
UVFELs and high average power FELs, the former due to its short wavelength and 
therefore a small diffraction rate, the latter due to extremely high (= 100 MW) 
intracavity average power. This problem could be solved by designing a long optical 
cavity with Z,) = 0.3, allowing space for diffraction to make the beam spot on the 
mirrors large. Another way to make a large spot is to use a cavity with a short 
(Z) < 0.3) Rayleigh length. When the FEL has limited space for the optical cavity (on a 


ship, for example), the second method of mirror protection is particularly attractive. 


As an example of the effect of the Rayleigh length on mirror power density, 
consider a 1 MW FEL at A=1 um and Q =5. Simple calculations based on (60)-(63) 
show that the mirror power density in a 20 m long cavity with Z)=0.2 is about 
1500 kW/cm?, which is quite large. However, if we reduce the Rayleigh length to 
Z, = 0.02 and shorten the cavity length to 10 m, the mirror power density is reduced to 
about 600 kW/cm®*. That is, we can reduce the power density on the mirrors by more 
than a factor of two and cut the cavity length in half by reducing the Rayleigh length by 


a factor of ten. 


Reducing the Rayleigh length is accomplished by making the mirrors larger and 
their radii of curvature smaller (recall Z,) = (2S5R - S*)/2L), allowing the light to cover 
a large mirror while focussed to a small waist at the center of the cavity. The stability 
cntena for a symmetric optical cavity, O < (1 - S/A)(1 - S/R) <1, reduces to ie 20 
(Ref. 7). Therefore, as long as there is a Rayleigh length (i.e., 2) > 0), the optical 


Cavity will be stable. 


In section G of Chapter Il, we discussed the physics of the FEL interaction in 


the Gaussian optical mode. We discussed how the optical field amplitude lag! is 
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greatest at t=7,,, which causes the phase space separatrix height to increase as 
t:0— 1, then decrease as t: 1, > 1. This effect, which is proportional to w~%(t), 
complicates the phase space paths and will be exacerbated by a small Rayleigh 
length. The short Rayleigh length implies a small mode waist where the optical field 
amplitude will be much larger than at the ends of the undulator. We also saw how the 
phase of the Gaussian mode caused the separatrix to shift to the right as t: 0 — 1. 
This requires the initial phase velocity for maximum gain to be increased from vo = 2.6, 
the value found when using the plane wave approximation. For Z,) = 0.1, the phase 
shift for on-axis electrons is Ago = o(t=1) — o(t=0) = 0.92, and is smaller for off-axis 
electrons. With a small Rayleigh length, the wavefront curvature has to be larger in 
order that the light be focussed to a small waist and rapidly expand to cover the 


mirror. 


Operating the FEL with a short Rayleigh length is not the optimum design. The 
previous paragraph showed how the interaction of the electrons and the light changes 
with Rayleigh length. This concept has never been studied before now, precisely 
because the design is not optimum. in this chapter we will show the consequences of 
using a short Rayleigh length. Since we will analyze the case Z) <1, we cannot 
expand (60) in 1/(Z,)) as we did in Chapter II, but must use simulation instead. We will 
analyze the UVFEL of Chapter IV first, both the single pass gain G(t) and the gain 
evolution G(n) over several passes when in steady state. Next we present a MW- 
class FEL design based on the short Rayleigh length concept and simulations of 


same. 


B. UVFEL WITH SHORT RAYLEIGH LENGTH 


We first test the short Rayleigh length FEL concept by simulating the LPC 
UVFEL, which has parameters listed in Table 2. 
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lis Single Pass Gain 


We study the single pass gain G(t) in the UVFEL using transverse multimode 
simulations like those discussed in the previous chapter. In this case, however, since 
only one pass is simulated, cavity mirrors do not enter. The optical field is initialized 
as a Gaussian mode and can only evolve by mode distortion caused by the electron 
beam within the single pass. Contrast this with the multiple pass simulations 
presented in the previous chapter in which the optical field, entrained by the mirrors, 


evolves self-consistently over many passes. 


Plotted in Figure 22 is the single pass gain evolution G(t) along the undulator 
for four different electron beam radii: (a) o, = 0.2, (b) o, = 0.45, (c) og = 0.9, and (d) 


Sng cb; 
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Figure 22. Gain along the undulator G(t) for electron 
beams with radius (a) o, = 0.2, (b) o, = 0.45, (Cc) 
6, = 0.9, and (d) a, = 1.6. Rayleigh length ts 2) = 0.3, 


mode waist radius iS Wo = 0.5. 
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The Rayleigh length is Z9 = 0.3 giving a mode waist radius of wy = 0.5 (see (63)). The 
expected LPC UVFEL electron beam emittance ¢ = 0.43 is again. Electron beams (a) 
and (b) have radii smaller than the mode waist radius (i.e., 6g < Wo) so that they are 
enveloped by the light for the duration of their interaction. The radii of electron beams 
(c) and (d) are greater than the mode waist radius. Because of this, their final gain 
G(t=1) is reduced compared with beams (a) and (b). Also note that because of the 
large angular spread oy, = 4.6 for electron beam (a), the final gain is smaller than that 
for beam (b). We saw in section D of the last chapter that electron beam (b) with 


Oe = 0.45 and od, = 0.9 was nearly optimum for Z, = 0.5 (optimum was o, = 0.4). 


The most encouraging aspect of the plots in Figure 22 is that although the final 
gain differs among the four electron beams, all show nearly exponential gain (see 
(58)) throughout the length of the undulator. This point is better illuminated by Figure 
23. In this figure we normalized the curves G(t) from Figure 22 to their respective final 
gains G(t=1) in order to compare the structure of the curves more closely. Figure 23 
shows convincingly that all four of the electron beams simulated resulted in similar 
gain evolution throughout the length of the undulator, even those ((c) and (d)) with 
Ge > Wo. We see that the gain evolution starts in all four beams after a small delay of 


t = 0.2, which is just the time it takes for the electron bunch to form (Ref. 3). 


In Figure 24 we have plotted G(t) for the same four electron beam radii, but 
with the Rayleigh length reduced to Zz) =0.1. This means Wy = 0.3 and only electron 
beam (a) with o, = 0.2 fits inside the optical mode throughout their interaction. For 
this reason, beam (a) has the highest final gain. Electron beam (b) has final gain 
reduced by nearly a factor of two over that in Figure 22 because it no longer fits inside 
the optical mode at the waist. Beams (c) and (d) also have final gain reduced, but 
only beam (d) has final gain less than 10%. This figure shows that in a single pass, 
where the electron beam cannot couple to higher order modes, it is more important to 


have the electron beam enveloped by the light than to have a small angular spread. 


64 


G(1)/G(t=1) 





Figure 23. Normalized gain along the undulator 
G(t)/G(t=1) for electron beams with radii (a) o, = 0.2, (b) 
6, = 0.45, (c) 6, =0.9, and (d) o, =1.6. Rayleigh 


length and mode waist radius same as in Figure 22. 


For Figure 25, the Rayleigh length was further shortened to Zp) = 0.02. In this 
Case Wo = 0.14 and none of the four electron beams fit inside the optical mode at the 
waist. The final gain for electron beams (a) and (b) is less than half their value at 
Zo = 0.1 in Figure 24. However, there is still gain at this very short Rayleigh length, 
provided the electron beam is small enough. We see that oc, = 0.9 or 1.6, Curves (c) 


and (d), are clearly too large. The final gain G(t) > 20% for the two smailest electron 


beams. 
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Figure 24. Gain along the undulator G(t) for electron 
beams with radius (a) o, =0.2, (b) o, = 0.45, (c) 
0, = 0.9, and (d) 6, = 1.6. Rayleigh length is Z, = 0.1, 


mode waist radius iS Wo = 0.3. 


In order to better compare the evolution of gain G(t) in Figure 25, the curves 
are normalized to their respective final values, as in Figure 23. The result is plotted in 
Figure 26. The effect of the short Rayleigh length is clear: single pass gain increases 
greatest where the fields !al are strongest, at the beam waist t= 0.5. The short 
Rayleigh length causes the fields at t= 0.5 to be much Stronger than at any other 
position. This can be seen in electron beams (a), with o, = 0.2, and (d), with 
6, = 1.6. Electron beam (a) is focussed at t = 0.5; it was designed to take advantage 


of the field strength at the center by concentrating electron beam current there. 
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Electron beam (d) has a radius eight times that of beam (a). It is so large that gain is 
poor throughout the interaction at this Rayleigh length. What little final gain exists, 
however, was enhanced at t = 0.5 where optical fields are strongest. Indeed, we see 
in Figure 26 that there was net absorption for beam (d) until the field strength started 
to increase at t= 0.4. Also evident in the figure is that electron beams (b) and (c) 
have radii which allow gain to be distributed more uniformly. That is, they are larger 
than beam (a) so gain does not increase as much in the center of the undulator, but 


are smaller than beam (d) so that gain increases more at the ends of the undulator. 
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Figure 25. Gain along the undulator G(t) for electron 
beams with radius (a) o, =0.2, (b) o, = 0.45, (Cc) 
6, = 0.9, and (d) 6, = 1.6. Rayleigh length its Z) = 0.02. 


mode waist radius tS Wo = 0.14. 
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Figure 26. Normalized gain along the undulator 
G(t)/G(t = 1) for electron beams with radius (a) o, = 0.2, 
(b) o, = 0.45, (c) co, = 0.9, and (d) o, = 1.6. Rayleigh 


length and mode waist radius same as in Figure 25. 


To summarize, we have found that using a ven short Rayleigh length degrades 
However, gain is still significant for the LPC UVFEL at a Rayleigh length as 
short aS Z) = 0.02. For o, = 0.2, the final gain is G = 37%, while for o, = 0.45, the 
final gain is G = 23%. The gain is only senously degraded for electron beams with 
G,. = 0.9 and o, = 1.6, which have radi two and four times greater than optimum, 
respectively. We also found that the most important region for the interaction is the 
center of the optical mode where the fields are strongest, even if the electron beam is 


much larger than the mode waist. This makes sense because it ts the field amplitude 
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lal which drives bunching in the pendulum equation (36). 


When this research was started, however, it was thought that the interaction 
may "turn off" at the center of the undulator, meaning that section of undulator could 
be replaced by a drift space, making an optical klystron, which would improve gain. 
The optical klystron consists of two undulator sections with a drift space at the center 
(Ref. 3). The first section, called the modulator, starts the electron bunching process. 
The drift space lets the electron momenta continue bunching in the absense of the 
undulator field. The second undulator section, called the radiator, causes the bunched 
electrons to radiate coherently. Instead of a drift space, the same effect can be 
obtained in a shorter distance with a dispersive section, which is just a magnet that 
bends the electron beam out of and back into the beam path. Slower electrons bend 
through a smaller radius than the more energetic electrons. Since the faster electrons 
travel a greater distance through the dispersive section than the slower electrons, they 
are bunched upon exiting. It was found, however, that we do not need to resort to 


exotic klystron undulator designs to get the desired gain. 


2. Steady State Gain 


Having studied how the gain G(t) evolves within a single pass, we will now 
simulate several passes and follow the gain evolution G(n). Recall that G(n) follows 
the final gain G at t=1 after each pass n. We use the transverse multimode, 
multiple pass simulations used in the previous chapter (see, for example, Figure 20) 
where the optical mode evolves self-consistently in the optical cavity over many 


passes. 


The electron beams (a), (b), and (c) from the previous figures were simulated. 
The simulations start from weak fields a) = 0.01 and evolve until the limit of weak 
fields lal =n is reached. At this point the final gain is in steady-state, by which we 
mean that single pass gain G(t=1) has nearly the same value over several passes 


and G(n) is nearly constant. It is this steady-state value of the gain evolution G(n) 
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that is plotted versus the Rayleigh length Zp in Figure 27. 
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Figure 27. Steady-state gain G versus Rayleigh length 
for electron beams with radius (a) 6, = 0.2, (b) o, = 0.45, 
and (c) 6, =0.9. We find G2 60% at Z)=0.01 for 


beams (a) and (b). 


In Figure 27, we see that the nearly optimum electron beam, (b), has the 
highest steady-state gain at all Rayleigh lengths explored. This is in contrast with the 
single pass studies in which electron beam (a) had highest final gain for Zp, < 0.3. 
This electron beam is evidently more efficient in coupling to higher order modes than 
any of the other beams explored. Perhaps this is because of the better balance 
between it’s beam radius o, = 0.45 and angular spread o, = 0.9 in their contribution to 
the electron phase velocity spread (i.e., it is better “matched” (Ref. 3)). The maximum 


gain for all three of the beams occurs at Rayleigh length Zp = 0.1. The most important 
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feature to note in the figure shows the goal of this chapter: steady-state gain is 
substantial even at Z) = 0.01. For electron beams (a) and (b) we have G(n) 2 60%, 


while that for beam (c), twice the optimum radius, is still a respectable 12%. 


Our conclusion is that the FEL with a short Rayleigh length may be able to 
operate above threshold with respectable gain while spreading the optical mode at the 
mirrors. The trade-off is slightly lower gain for a much shorter optical cavity with 


reduced power loading at the mirrors. 


C. MW-CLASS FEL WITH SHORT RAYLEIGH LENGTH 


Encouraged by success in modelling the UVFEL with a short Rayleigh length, 


the concept was used to design a MW-class FEL. 


1. MW-Class FEL Design 


Our goal in designing a MW-class FEL is to derive a set of parameters which 
will provide sufficient gain within an optical cavity of length 10 m. This length was 


chosen because it fits within the beam of a modern naval cruiser or destroyer. 


FEL gain, discussed in section E of Chapter Il, is dependent upon the 
dimensionless current density / (see (57) and (58)). Maximizing gain therefore 
requires that we maximize /. Recall that the definition of the current density in the 


linear undulator is 
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where & = K*/(2(1 + K*)). Since N = L/A, we can re-write (69) as 
2 | K*(Jo(&) - J4())*L° 
| 70 
; me? Ao ~ BS, 


where we have separated the parts of / according to their source. Following the 


fundamental constants is the undulator contribution K*(SI\SL AG where we have 


71 


shortened the Bessel function term in (70) to JJ. Since K = 1 in an FEL, which also 
fixes the Bessel terms, we are left with L and A, as designable undulator parameters 


(Ref. 3). The electron beam contribution is p/y°. The charge density p depends on the 


peak current / through 


Q 
Se Til 
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where Q = / x (/,/c) is the charge per micropulse (not the cavity quality factor), r. Is 
the electon beam radius, and /, is the length of the micropulse. The electron beam 


energy ymc* determines y. Therefore we can write 


pe eeaine 
oc ——-————__, (72) 
Ao ie Y 
Recall that the emission wavelength is given by 
2 
Na ena (73) 


ai 
With K =1 this becomes A = A)/7*. This means that if we fix the wavelength A, then 
Xo and y are not independent. Therefore, the remaining free parameters are 


ite anions 


Our philosophy in choosing a design was to use a short undulator and use 
peak current to maximize /. The reasons for this are twofold. First, we are using an 
optical cavity 10 m long with a short Rayleigh length. The undulator has to be short in 
order to fit around the optical mode at the center yet not “scrape” any of the light, 
which would damage the magnets. Second, a short undulator means the FEL will be 
less susceptible to electron beam quality. This is because the gain spectrum 
bandwidth, which is proportional to 1/(2N), increases for decreasing undulator length 
L = Niy. This was discussed in section E of Chapter Il and is important to the design, 
because as we demand more current from the accelerator we expect the beam quality 


to suffer. 


fe 


Accordingly, we chose an undulator length L = 1m. The undulator wavelength 
was taken to be Ag =4 cm, which means N =25 periods. Since the wavelength 
X= 1 um is fixed by the ship self-defense application, the electron beam energy is 
ymc* = \Ag/Amc? = 100 MeV. For the remaining parameters we chose / = 600 A, 


fe =0.3 mm, and /,/c =3 ps. This results in a total charge per micropulse 
Q 


i x (l/c) = 1.8 nC. If every RF bucket of an f = 500 MHz accelerating field is 
filled by the electron injector, the duty factor is (/,/c) x f =0.0015. The average 
current is then 7 =0.0015x/=0.9 A. This results in an average electron beam 
power of P, =7V=0.9A x 100x10®° V = 90 MW. An efficiency of 1.1% will be 
needed to extract 1 MW from this beam. These design parameters resulted in j = 12. 


All parameters are summarized in Table 3. 


For electron beam quality we anticipate a trade-off between high peak current 
and electron beam quality. We assumed an energy spread Ay/y=0.5% and a 
normalized emittance €, = 20m mm-mrad. These choices resulted in a = 1.5 and 
€=0.27. The dimensionless electron beam radius is o, =0.5 so that the angular 
spread must be o, = 0.3. Noting that N is included in the numerator of (18) and (56), 
we see that the small values for 6g, 6g, and € reflect our choice of a short undulator. 


The beam quality parameters are also summarized in Table 3. 


2. Simulations of the MW-Class FEL 


Now that we have derived our parameters for the MW-class FEL, we use 
simulations to make sure it works. Transverse multimode, multiple pass simulations 
using the parameters in Table 3 were used to obtain the steady-state gain at different 
Rayleigh lengths. An example of these simulations at 2) = 0.1 is shown in Figure 28. 
Steady state gain is G(n) = 26%, plotted at lower right, and power P(n) ts increasing 
approximately exponentially. The optical field evolves over 20 passes yet s still weak 


(la! < x). Bunching is not evident in the phase space plot at lower center. 
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Beam Energy 

Peak Current 

Average Current (@500 MHz) 
Micropulse Length 

Charge per Micropulse 
Electron Beam Radius 

EN 

Ay!y 

Undulator Type 

K 


N 
ho 
r 


100 MeV 
600 A 
0.9 A 
3.0 ps 
ee 

0.3 mm 

20x mm-mrad 
0.5% 
Linear 
1:0 
25 





Table 3. Design parameters for a MW-class FEL suitable 


for ship self defense. 
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Figure 28. Transverse multimode simulation of MW-class 


FEL. Gain evolution G(n) is in steady state. 


The electron positions overlaid on the plot of the optical field amplitude between the 
mirrors |a(x,t)| show that the electron beam does not expand nor focus as 7:0 —> 1. 


This is because the angular spread o, = 0.3 is relatively small. 


Many simulations such as that in Figure 28 were used to find the effect of the 
Rayleigh length on steady-state gain. The result is plotted in Figure 29. As in our 
previous analysis of the UVFEL, maximum gain G = 23% for this FEL occurs at 
Z, = 0.1 ( see Figure 28). Even at the short Rayleigh length 2, = 0.01, the steady- 


state gain G = 18% is sufficiently large that simulations indicate the design will work. 


aS 


02 


O15 


0.1 


0 0.05 O01 0.15 0.2 0.25 0.3 
Z0 


Figure 29. Steady-state gain G versus Rayleigh length 
Zo for MW-class FEL. Peak gain G = 23% occurs at 
Zo = 0.1. At 2) = 0.01, gain is G = 18%. 


CONCLUSIONS 


Using the LPC UVFEL parameters we showed that the concept of the short 
Rayleigh length FEL was viable. This had never been researched before. Based on 
those results, parameters for a MW-class FEL for ship self-defense were denved. This 
FEL design was simulated and steady-state gain was found to be G = 18% at 
It was presented at the workshop entitled “Navy MW-Class SSD FEL 
Concepts” at TJNAF in Newport News, Virginia, which was sponsored by the Navy 


High Energy Laser office at SPAWAR. It is one of only two designs that survived the 
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workshop. 


To improve on our design, the electron beam radius should be optimized. At 
O, = 0.5, we have a small 6, = 0.3 which suggests it would be advantageous to focus 
the electron beam further. Research should also determine what affect the Rayleigh 
length has on FEL efficiency. At short Rayleigh lengths, many electrons are outside 
the optical mode at the waist and do not participate in the interaction. Efficiency might 
be degraded. This is the trade-off made in choosing a short Rayleigh length over 
more traditional FEL designs. As we noted above, we must extract 1.1% of the 
electron beam energy in order to obtain 1 MW average power. Methods of improving 


efficiency should also be studied. 


The design philosophy of using a shorter optical cavity and undulator requires 
high peak and average current. It shifts the technology burden from the mirrors to the 
accelerator, where improvements in high average current generation and acceleration 
must be made. Another area of accelerator physics that needs much further research 
is in electron beam transport. For our laser to work on a ship, we have determined 
that we must use an energy recovery architecture similar to that proposed for the LPC 
UVFEL in order to increase overall system efficiency, among other reasons (Ref. 9). 
This reduces the power consumption of the FEL by recirculating the electron beam 
through the accelerator, where it is decelerated, giving up its energy to the RF 
accelerator field. Recirculation requires the bending of a high current, high energy 
electron beam through two 180 degree bends, but, at the time of this writing, there is 


some question as to whether or not that is even possible. 





Vi. COMPTON BACKSCATTERING OF LASER LIGHT 


A. PRELIMINARIES 


We turn our attention now to the scattering of laser light from a high energy 
electron beam. The electron beam scatters high energy photons out of an incoming 
laser field. The process is similar to the FEL in that the incoming laser pulse can be 
thought of as an electromagnetic undulator, but with extremely small wavelength. The 
process has been described classically and with relativistic quantum mechanics, 
assuming a plane wave approximation for the laser field (Refs. 11, 12, 13). Recently, 
a correction to the plane wave approximation was attempted, using quantum 
electrodynamics (QED) (Ref. 14). Derlet et. al. treated the interaction of free electrons 
with rays of laser photons focused by a lens to the interaction region (Ref. 14). The 
calculation presented in this chapter is the first attempt to correct the plane wave 
approximation to the laser beam by treating it as a diffracting beam. That is, we 


assume the incoming laser beam is a fundamental Gaussian mode. 


In the next section, we calculate the vector potential for the Gaussian optical 
mode. We approximate the Gaussian mode for long Rayleigh length in a manner 
similar to the analysis in section G of Chapter Il. In section C, we calculate electron 
wavefunctions in the field represented by our vector potential. We use those 
wavefunctions in section D to find the distribution and rate of photon emission. This 
will be the first time this calculation has been done without the plane wave 


approximation to the laser field. 


The notation used in this chapter is that of Sakurai (Ref. 15). We represent a 


four-vector as, for example, b,, with yp = 1,2,3,4 where 
by, = (Dy,D2,D3,b4) = (B,ibo) (74) 


where b,, D>, and bs, are real, corresponding to the x-, y-, and z-components of the 


three-vector Bb, and b, = iby is imaginary. We use a summation convention so that 


oo 


the scalar product of two four vectors 5, and c, is 
3 > 
nen = P56; —- DACA =D 0a Dolo ‘ (75) 
= 
We will frequently shorten the notation for a scalar product from bic, to b-c for 


brevity. In our summation convention, Greek indices such as un, v, etc. run from 1 to 4 


and Roman indices like /, /, k, etc. run from 1 to 3. We define the coordinate vector 


as 


xX, = (X,ict) (76) 


where X = X;X + Xof + X32. The momentum four-vector is 


Py = (D,iE/c) , (77) 


where p and E are the particle momentum and total energy, respectively. The 


wavevector of a photon is 
k, =(K,iwic) , (78) 


where I|k| = 27/4 = w/c for a free photon. Variables are not made dimensionless as in 


the previous chapters. 


We use Heaviside-Lorentz units throughout this chapter. The fields and 
potentials in this system are related to those in Gaussian units by the factor 1/V4n 


(Ref. 15). We will take the electron charge e to be negative e = lel. 


A typical experiment, and the one we will use in order to obtain order-of- 
magnitude estimates to certain terms, is the proposed “laser synchrotron source” 
(LSS) at the Naval Research Laboratory (NRL) (Ref. 11). They have proposed using 
a 41 MeV electron beam with micropulse length /,/c = 1 ps. The incident laser is a 
tabletop-terawatt (7°) device with wavelength A = 1 um. The laser peak power is 10 
TW in a 2 ps pulse, which means the energy per pulse is 20 J. The Rayleigh length is 
Zo = 7.9 mm and the mode waist radius is Wp = 50 um. The energy of the emitted 


photons is expected to be 30 keV (x-rays). We can then calculate some typical values 
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for useful combinations of parameters: eA ,/mc* = 0.4, hOQ/mc* = 0.06, and 
ha/me* = 2x10-®. We define Ay as the amplitude of the vector potential, m is the 
electron rest-mass, c is the speed of light in vacuum, fi is Planck’s constant divided by 
2m, Q is the frequency of the emitted photon, and ow is the frequency of the incident 
laser. Note that the combination eA,/mc* is analogous to the undulator parameter K 


in the FEL. 


B. VECTOR POTENTIAL FOR LASER LIGHT 


The fundamental Gaussian mode for some freely propagating scalar field U 
can be written: 
x?+y? 


Z Uy a ae (79) 
w(z) 





where U, is the amplitude in the center of the beam waist (Ref. 7). The mode radius 


w°(z) = zh + 21] (80) 
20 


where Ww, is the radius of the mode waist. The mode waist radius and Rayleigh length 


is 


are related by 





(81) 


The optical phase 9 is 


6 = “en 2) + Pasa dl (82) 


Zo 2A(z2) 


where A(z) = 2+ 2¢/z is the wavefront radius of curvature. The origin z = 0 is 


located at the mode waist. 
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Since our goal is to obtain a correction to the plane wave approximation for our 
scattering problem, we will find the form of (79) when it nearly represents a plane 
wave, Or Z) > Z. For this case, it makes sense to expand (79) in the small parameter 
(1/Z9) in a manner similar to the WKBJ approximation (Ref. 16). That is, we will 
expand both the amplitude and phase of (79) in (1/Z)) and approximate them by 


keeping terms only to first order in (1/Z,). 


Upon expanding the amplitude we find that 











U 34] U 2 9 2 
0, (wey) =O Key) ee ee (83) 
w(Z) Wo Zo 2 20 
Likewise we find for the phase 
,)2 
(== 1] . (84) 
20 Zo 
Therefore the WKBJ approximation for the Gaussian mode is 
U = =I : pe era Zoi Zoe (85) 
Wo Zo 2 


To lowest order in (1/Z)) we see that the Gaussian mode is nearly a plane wave. The 
only differences between (85) and a plane wave are a phase shift from k > k —- 1/Z, 
and an amplitude correction which decreases the amplitude for electrons off-axis. To 
an electron in the Gaussian mode, the phase shift appears as a shift in the laser 
wavelength A — A/(1 — 1/kz.) = A(1 + 1/kZ_). Therefore we can anticipate that the 
presence of the Gaussian mode should result in a shift in wavelength of the emitted 


radiation as well. 


For our calculation, we will take (85) to be our vector potential. Since this 


potential is real we will use the real part of (85) and write 


yh 1 As 
a 


0 


A =Aoa 





| cos| (k - 11z9)2-wl (86) 
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where Aj, with Wo absorbed, is the on-axis magnitude of A and 4 is the polarization 
unit vector. If the laser light is linearly polarized in the x-direction, as we will assume 
in this chapter, then & =(1,0,0). We will only consider on-axis (i.e., x? + y? = 0) 


electrons, so that we can write (86) as a four-vector as 


A,,(x) = (A(x), iAo) = Aya,cos(K,x,) (87) 
where a, = (2,0) and 
Ce Coad (88) 
Zo C 


is the modified propagation vector. The vector 6, = (0,0,1/Z),0) represents the phase 
shift discussed above. We note that A,(x) satisfies the Lorentz condition 
0A,,/dx, = K,A, = 0 and, since in the gauge we have chosen Ap = 0, the transversality 


condition V-A = K-A = 0 is also satisfied (Ref. 15). 


C. DIRAC WAVEFUNCTIONS 


The wavefunctions for an electron in a plane wave electromagnetic field were 
found in 1935 by D. M. Volkov (Ref. 17). We will follow nearly the same technique, 


but using our vector potential (87) for the Gaussian mode. 


The Dirac equation for an electron in an electromagnetic field characterized by 


A,,(x) is given by 





02 Heya EO (89) 


e ax, Nhe nh 


where the 4x4 matnces y, are called the “gamma matrices” (Ref. 15). They are 


defined by the requirement that they satisfy 


(Yue Yo? = Yute + Wy = Ope (90) 


where { , } is called the anti-commutator, and 6,,.= 1 if » =v, 0 otherwise, is the 


83 


Kronecker delta (Ref. 15). In our notation, the gamma matrices are 


0 —lOx i O 
Yk = wa 0 | » ¥4 = kl a ’ (91) 
where / is the 2x2 identity matrix and the o, are the Pauli spin matrices given by 
0 1 0 -i ie 
0, = [' 7 02 = |" 7 03= [4 mI (92) 


For example, we write y3 using o3 as 


00-10 
000i 

(momo O| an S 
0-100 


We also note here that the gamma matrices are Hermitian y; = y, (Ref. 15). The 
wavefunction in (89) is a four component object called a Dirac spinor (Ref. 15). Thus 


(89) is actually four coupled, first-order, differential equations for w. 


We can obtain the quadratic form of (89) by multiplying on the left by the 


operator 
"OX, fce'# RP 


with the result 








da 2ie, a ie dA, e@ mc? 
ee —~-——|w=0 . 95 
Oxo xiemnCemeoX,, ey oz bP . (eo 
Substituting (87) for A, we find that 
dA, 
— =K,A,’ (96) 
OX, 
where 
A,’ = -Apasin(KgX_) . (97) 


We can then rewnte (95) as 


0 oO 2ie d ie e* mc? 
Fe eee A* - ——ly=0 98 
OX, OX, NC Oeeie pee) tec? ir ns a 





where we have used y-K = y,,K,, etc. 


Following Volkov, we seek a solution to (98) of the form 
Px 

w=e ” o(k-x) (99) 
where @(k-x) is a Dirac spinor. We subject our solution to the initial condition 
yvyo u'S)(B)exp(ip-x/h) as x > -c<- for an electron coming into the interaction region. 
We recognize u'S)(f) as the free electron spinor, s = 1,2 is the electron spin index 
and p is the electron momentum (Ref. 15). The initial condition means physically that 
prior to entering the interaction region the electron is free. Substituting (99) into (98), 


we find that o(k-x) must satisfy 


S pk oie) + |S Ap ~ SB crkyira) - 








cs =O (100) 


where &=k-x, 6’ = 00/0&, and we have used k*=0. This first-order differential 


equation has the solution 


o(E) = exp 2. j Uk ee ai 





A2 
AD 8 A ea 101 
Apt Magia (on 
where u is a Dirac spinor. Invoking the initial condition, we find that ¢ — u as § > —- 
so that y — uexp(ip-x/h). Therefore, we will take u to be the free electron spinor 
u = Vmc*/EVu'*)(p) where V is the normalization volume and E is the total energy 


(Ref. 15). The wavefunction for the incoming electron is 


| . g 

@e (y:K)(y-A’) a ie 

26 J ok ee hee 
“x 


mc? : (s) : n 
s 
x EV u’’'(B)e 


Ap e@ At 
p-k 2c p-k 





jas (102) 


2) 


—O 


85 


We can show that this wavefunction satisfies the Dirac equation (89), as it 
must, in addition to being a solution to (95) (Ref. 12). Let D represent the Dirac 
equation operator such that the Dirac equation is Dy =0. Then DDy = 0 represents 
(95). (Recall that D was defined in (94).) Let y=Dwy,, where w,, is given by (102). 
Now since w;, reduces to the free electron spinor, which satisfies the Dirac equation, 
as X > -co, we know that y > 0 as x ~-~. Since DDw,, = Dy =0 and y > 0 as 
X99 — -co, we must have x = 0 identically. Therefore, Dw,;, =O and y,, satisfies the 


Dirac equation (89). 


To find the wavefunction for the outgoing electron, we must solve (95) subject 


to w > u'*)(B’)exp(ip’:x/h) as Xp > e. The result is 


Z _er (VK )YA) ger _ jet Ap _@ A® der’ 103 
VW out = exp 20! p'-k é ne! p'-k 2C p k C ( ) 


me Ya iP x 
me” (Ss) (R’ n 
u D e 

a ae 


x 





where the outgoing electron is specified by spin index s’, momentum p’, and energy 


ce 


We now turn to the task of evaluating the integrals buried inside the 
wavefunctions (102) and (103). The vector potential in the integrals of (102) and (103) 
is a function of K-x = k-x —8-x while the variable of integration is k-x. Let 


y = K-x = E-6-x and perform a change of variables to xy. We find 


0 1 1 
1+ —(8-x)|/dy = |1+ — + O|— 


dz = 
be) Zo 





ox (104) 








and we approximate the differential d& with the first two terms. Making the change of 


variables in (102) we wnte 
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x+d°X 
z 1 _||_e Veo ants) ee 
ao 1 ae 
Win oe || a J De J ee dy (105) 
_ xto-x , Clas 2 te eS 
+ € Ait cp ae” Ante dy’||| me u(S\B)e h 
fic 2. p-k 2c pk EV 














The correction term (1+ 1/kZ)) which we discussed in the last section appears 
explicitly in the electron wavefunction. Substituting the vector potential (87) into (105) 


and performing the integrations we arrive at 








Win = eX h + a 2 (vide) cos(k-x) (106) 
je|, |apl,., Ad} 1 || k-x . sin(2k-x) 
+ = Ao 22 sm x) ye || 5 + =e } 











_p:x 


Ya 
mc? (s) : h 
x | ne u”'(py)e 





where we have dropped divergent sine terms in the phase of the wavefunction. The 
divergence of these terms is not physical, arising because our vector potential 
describes a laser beam undamped at -—<- while we require that the fields associated 
with A,, vanish at -cc. Likewise, we perform the change of variables (104) and the 


integration in (103) with the result 
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ws 1 e || (:K)(a) 
= Becton ; 
Wout = eX 1 elle avy cos(k-x) (107) 
. 2 
ie ole CAg| 1 || k-x — sin(2k-x) 
anal A tony EES jp Ne 
Z fc | ° Bi x) 2C nal 2 i 4 
oe jP-X 
«| RE) um pre’ ? 


for the outgoing electron. 


We now have the wavefunctions for the electron entering the interaction region 
(106) and for the electron leaving the interaction region (107). The difference between 
these wavefunctions and the Volkov wavefunctions is in the factor (1 + 1/(kZ9)) in the 
exponential and the presence of the modified propagation vector K,, instead of k,, 
(Refs. 12, 13, 17). As we expected, by choosing a Gaussian mode which differed only 
slightly from the plane wave (i.e€., Z) > Zz), our wavefunctions differ only slightly from 


the Volkov wavefunctions. We now use these wavefunctions to find the effect of the 


Gaussian mode on photon emission. 


D, PHOTON EMISSION 


For the calculation of the differential transition rate in this section we will use 
“natural units” N= c= 1. This aids the calculation in that we will not have to keep 
track of all of these fundamental constants. We will wnte our results at the end of the 


section with these constants included so that real numbers may be calculated. 


The amplitude for an electron transition from y,, tO Wo, with the emission of a 


photon 1's 
M = -{j,Ald*x (108) 
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where the integral is over all space and time (Ref. 15). The interaction vector potential 


A") for the emission of a single photon is 


Aine ala 


m= NV payne” (109) 


where /,, is the emitted photon propagation vector and «!"), r = 1,2, is its polarization 
vector (Ref. 15). These are given by 

ee) (110) 
where Q is the frequency of the emitted photon and 


ae = (NG (6) (111) 


is chosen such that (a) e@) 7717) form a right-handed coordinate system in three 
dimensional space (Ref. 15). The interaction of the electrons with this unquantized 
vector potential (109) for photon emission is rigorously equivalent to their interaction 


with the quantized electromagnetic field (Ref. 15). The transition current j,, is 


Julx) = 16 Vout Yuin (112) 


where y, is given by (91) and y=w"y, is the adjoint spinor (Ref. 15). Substituting 


(109) into (108) we can write the transition amplitude as 
i ~i(I-x) 44 
M= saz Jintrele Ht va, x . (113) 


To find the transition current j,,(x) we substitute the wavefunctions (106) and 


the adjoint of (107) (recall M = c = 1) into (112) with the result 
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jul) = len] Zea] A elle Ps (114) 


COsk:x/¥,, 


eA ‘ay 
x4 O'S )(B)exp| —"]1 + all vem ass 
2 p':k 





yKy-a 


ai cosk-x| u'S)(B) 


eA 1 
x exp yr 1+ in 
0 











x exp|iasink x — IB(2k-x + sin2k-x)] 


where we have defined 











1 ap asp 
= €@A,|1 + ——!|—~— - 113 
a= CAg|1 + 2s ae (115) 
and 
e-Ag 
poe eee ee (116) 
8 kZ2)|| pK pk 





We will now simplify and make some approximations to (114). Using the 


property of gamma matrices that 
(y-b)(y-b) = b* (117) 
which is a special case of 
(yO)(y-c) = b-c - (y-C)lyb) (118) 


we can expand the real exponentials in (114). The expansion yields 
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PAo|, _ 1 | (ya)(yK) ono) _ 
oo : ir =n SER cosk:-z| = (119) 


eA 
ieee 





ee sa Yeas) cosk-x 
kZpo p -k 


where we have ignored terms of higher-order in (1/Z 9) and those proportional to 
(CAg/(2p’ )*(1/kZ_) = 10°'° and higher. We used the NRL parameters discussed in 
section A to calculate these estimates. We can simplify (114) by writing the complex 


exponential involving a and B in terms of Bessel functions using the generating 


function 


e/xsine =e y eine adoG) (120) 


M=—co 


where m is an integer (Ref. 18). We find 


expliasink-x - iBsin2k-x]}= YF ell™2Mkxy ey (B) (121) 


{7} =—co J? =—00 


= > eh yy Soa oe) 
n=—oo 


q =-9©o 


= e/@%* Go(q,a,B) 
q=-<0 


where G,(q.0.B) = > Yaqsn)(@)YUp(B). This term (121) leads ‘to emission into 
Nh=~—oc 


harmonics where q is the harmonic number. We then incorporate these simplifications 


and re-write (114) as 


j(x) = ele Ve T Golq.c.,Bye'OF *2(9-PikY* (122) 
qQ>-— 


eA enisgs 


x op) ice ek 


CA 1 
Ae 





Te cosk . u'S)(P) 
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Substituting (122) into (113) to get the transition amplitude we find that we can 


write 
M=Mo)+M,+Mb, (123) 
where 
Pe oto 4{ £4 4y @-ill-p+p-2(q-B)k)-x 
9 = —le s0V V Ev 7 Le ae (124) 
a "(pve u'S\(B) Go(q,0,8) 
and 
=-— 4 i(/-p+p’ -2(q-B)k)-x 
= -ie9/ mWVS EV EV, ave x Jo _ "ee 
eA 1 S aveK vel) vk yg 
2914, lS py) YEAVRY i VoL (s) 





x Go(q,a,B)cosk-x 


and M, « (eA,/m)*(w/m) = 10-’ can be neglected in comparison with My. For the 
transition probability, we need to find IMI* = IMol?+MgM,+M Mi +1M,|*. The 
term IM,|* « (eAg/m)*(w/m)* = 107'S and can be safely ignored in comparison with 
IM,l*. Also, the cross-terms like M{M, « (eAp/m)(w/m) = 10°, so that ignoring these 
terms when calculating |M|* will introduce an error of one part in 10°. Therefore we 


will use |M|* = IM)|*. Carrying out the trivial integration in (124) we find 


IM.12 = e2(2n)8—_ __™_ 15t4y_ p49 -21g-B8)k)I2 126 
M, Cxcn) === nT EV ey (/-p+p’ -2(q-B)k) (126) 


x lop ye u's)(p)l? G8 (q.0,B) 


for the specific harmonic number gq. We note that in order for the energy-momentum 


conservation relationship inside the 6-function to be true, q must be a positive integer. 


oe 


The total transition probability is then > IM, IF. We write the square of the delta 
q=1 


function as 


VT 
(2n)* 





18)(I-p +p’ -2(q-B)k) |? = 8(I-p +p’ -(q-28)k) (127) 


where V and 7 are the interaction volume and time (Ref. 15). 


The transition rate per unit time |Mq |*/T into a group of final electron and 


photon states is 


2 
_ IMG! 
i 





dW, (128) 





Va3 |p"! 
(2n)° 





Vas 71 
(2n)° 


where the two differential factors are the density of final states for the outgoing 
electron and emitted photon and integration over p’ and 7 is implied (Ref. 19). 


Substituting (126) into (128) we find 
2 
dw, = + © ™ s4\j_p 4p -(q—28)k) lo '"(p’)y-e usp) 2 (129) 
x Gé (q,a,B) dlp"! aS I7 


and we see that the dependence on V and T disappears, as it should. Performing the 


integrations we are left with 








-1 
dw, ee? Qmm| a Cay ac 2 
peered ome ao 4 Es 2 lo ly lS) py) 130 
a (ome 2 EE’ a0. ‘es: =) (B’)y (2) (130) 
x Gé (q.a.B) 


which ts the differential transition rate for a photon emitted into the solid angle element 
dx. The derivative term in brackets is the result of the [7) = Q integration over a delta 
function whose argument is itself a function of the photon emission frequency (Ref. 


19). We note that (130) is subject to the energy-momentum conservation relations 
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Pu = Py + 2(q = ly (131) 


which implies 


Pp’ =p +2(q-B)k-T (132) 

P= £424 —-pjo— oO (133) 
and because 

F = (ipl? +m)? , (134) 


E’ itself depends on Q through (132). Equation (131) is Lorentz invariant; it holds in 


all inertial reference frames. 


The remainder of the calculation will be easiest if we use the "laboratory frame" 
defined, for historical reasons, as that in which the incoming electron is at rest. This 


implies that BP = 0 and p,, = (0, im). Thus (132) becomes 
Pp’ = 2(q - B)K -T (135) 
and the outgoing electron energy is, substituting (135) into (134), 
Ya 
E = la(q — B)*w* + 2? - 4(q - B)wQcosé + m?| (136) 


where 6 Is the angle between k and 7. Since we know K = (2n/A)2, 6 measures the 


angle of photon emission from the z-axis, as shown in Figure 30. 


We can find the frequency 2 of the emitted photon by squaring both sides of 


(132) to obtain the invariant equation 
2(q - B)p:k - I-p - 2(q - B)l-k =0 (137) 


where we have used p* - p’* = 0. Computing the scalar products in (137) in the lab 


frame and solving for 2 we find 


Oe 2(q - B)w | (138) 
Teo By (1 ~ cos6) 


This frequency is corrected in that it is slightly decreased from the plane wave 
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approximation, since B « Aé(1 + 1/(KZ9)). The frequency is shifted by the factor 
(1 + 1/kZ9), which is the shift in frequency of the incident laser light as seen by the 


incoming electron, as discussed in section A. 





Figure 30. Coordinate system showing kinematics of 


Compton scattering (after Ref. 15). 


To compute a, (115), and 8, (116), we compute the scalar products: 


p-k=-mw , (139) 
pk =p’ k - Ew = of2(qg - B)w — Qcos6] , (140) 
ap=0 , (141) 
and 
a:p’ = &-p’ = -wsinOcoso . (142) 


Substituting these into (115) and (116) we arnve at 


( | 
ce eAo|1 - 1 SINBCOSO 


— | (143) 
kZ_ | 0(1 - cos8) - m 


where o 1s the angle between the projection of T onto the x-y plane and the x-axis (see 
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Figure 30) and 


e-A- 


8am 


1 


ed lice, Q(1 — cos) 
kZ 


m — Q(1 — cosé) 





B = (144) 








and we see the dependence of B on the emitted photon frequency 2. 


We are now in position to calculate the derivative in (130). Using (136), (138), 


(144) and a lot of algebra we find that 


0 ; 2qam 
— (2 = 
50 (Q+ E’ + 2Ba) FO (145) 





We recognize this term, apart from the harmonic number q, from ordinary Compton 


scattering (see Ref. 15, page 228). 


We can now write (130) in the lab frame as 


dWa QO 


oO 
aD An 
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qa 





lope yy-eu'S)(p) I? G2 (q,a,8) (146) 








where a = 1/137 is the fine structure constant (Ref. 15). The last thing we need to do 
is calculate the term luy-eul*. Since the spin of the initial and final electrons is 
unobserved, we average over the initial spins and sum over the final spins with the 
result 


E’ -m 
2m 





> ¥ lop yy-eulS\(pyl2 = (147) 
SS 


where we have used "Casimir’s trick" in order to compute the spin sum using the trace 


method (refer to pages 191-192 in Ref. 15). 


Our final result, for the emission rate of a photon of frequency 2 into the solid 


angie element df around 8, including /f and c, is 


dW, Q 


dy 4n 


2 
gw 














’ 2 
Emel gs (q.0,B) . (148) 


Recall that for the entire spectrum, (148) must be summed over the harmonic number 
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q. We summarize our other results, including f and c, as well: 


0- Se ae = eo (149) 
1+2(q= pa" — cosé) 


Va 
F = lag — B)*Fw* + HPQ? — 4FF(q - B)wQcosé + mec4| (150) 
1 sin6coso 
= €Ay|1 + — | ——— Tot 
i : kZ9 | AQ(1 — cos@) — mc? a 








er Ad |, a @) 
———|'{ + ———|(1 — cos 
seas kZo 
(152) 


te zs" a cos6) 


N 
On 





and we have used (138) to write B in terms of the scattering angle and laser frequency 
w. Our Gaussian mode correction (1 + 1/kZ,) is evident in each of these equations 
either explicitly or through 8B. We see that both the angular distribution (through Go in 


(148)) and frequency (149) can be altered by merely focussing the laser beam. 


If we specialize now to backscattering, 8=n, we see that n — 0 and the 


Bessel function sum collapses: 
Go(q 1=0,8) YY 2(n+q)(0)Yp (B) = (-1)7 J, (B) (153) 


since only Jp(0) > 0 implies n = -q and we have used the identity (Ref. 20) 
Jin (X) = (-1)" dp (x) (154) 
Therefore we find that the transition rate becomes 
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qu 


OW 
cma 


dx 4n 
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The argument of the Bessel function becomes 


on, 











0 
1+ 1/kz 
15 © 0 
a (156) 
moo + ~ 1 + 11kz9 
2mc 
and the emitted wavelength is 
oe “ig — ple (157) 
AP ACS) Saeee 
mc 


Since (fw/mc*) = 10~°, we can approximate Q by the numerator in (157). Writing this 


in terms of the wavelength we find 


vege 
N= oop = gl + Bia) (158) 


since 8B = 0.04. We see that the wavelength of the emitted photon is shifted by 


(1 + 1/kZ_) by the Gaussian mode, as expected. 


= CONCLUSIONS 


In 1935 D. M. Volkov derived the Dirac wavefunctions for an electron in an 
electromagnetic plane wave. We made a correction to those wavefunctions by solving 
the Dirac equation in the external field of a laser with a realistic, finite Rayleigh length 
Zp. We then used those wavefunctions to find the rate and angular distribution of 
photon emission in this field. We found that both the angular distnbution and 
frequency of the emitted photons can be altered by changing the Rayleigh length, that 


is, focussing the laser. 


Further research should add the finite length of the laser pulse to the 
calculation. This will give the emission lines a width in frequency corresponding, 


through the Founer transform, to the pulse length. 
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